Study of the A(e,e'7r+) Reaction on X H, 2 H, 12 C, 27 Al, 63 Cu and 197 Au 

X. Qian, 1 T. Horn, 2 ' 3 B. Clasie, 4 J. Arlington, 5 R. Asaturyan, 6 F. Benmokhtar, 2 W. Boeglin, 7 P. Bosted, 3 
A. Bruell, 3 M. E. Christy, 8 E. Chudakov, 3 M. M. Dalton, 9 A. Daniel, 10 D. Day, 11 D. Dutta, 12 - 1 L. El Fassi, 5 R. Ent, 3 
H. C. Fcnker, 3 J. Ferrer, 13 N. Fomin, 11 H. Gao, 1 K. Garrow, 14 D. Gaskell, 3 C. Gray, 9 G. M. Huber, 15 M. K. Jones, 3 
N. Kalantarians, 10 C. E. Kcppcl, 8 ' 3 K. Kramer, 1 Y. Li, 10 Y. Liang, 16 A. F. Lung, 3 S. Malace, 8 P. Markowitz, 7 
A. Matsumura, 17 D. G. Meekins, 3 T. Mertens, 18 T. Miyoshi, 10 H. Mkrtchyan, 19 R. Monson, 20 T. Navasardyan, 19 
G. Niculescu, 13 I. Niculcscu, 13 Y. Okayasu, 17 A. K. Oppcr, 16 C. Pcrdrisat, 21 V. Punjabi, 22 A. W. Rauf, 23 
V. M. Rodriqucz, 10 D. Rohe, 18 J. Sccly, 4 E. Segbcfia, 8 G. R. Smith, 3 M. Sumihama, 17 V. Tadevosyan, 19 
L. Tang, 8 - 3 A. Villano, 24 W. F. Vulcan, 3 F. R. Wesselmann, 22 S. A. Wood, 3 L. Yuan, 8 and X. Zheng 5 

1 Triangle Universities Nuclear Laboratory, Duke University, Durham, NC, USA 
2 University of Maryland, College Park, MD, USA 
s Physics Division, TJNAF, Newport News, VA, USA 
^Laboratory for Nuclear Science, Massachusetts Institute of Technology, Cambridge, MA, USA 
5 Argonne National Laboratory, Argonne, IL, USA 

6 Yerevan Physics Institute, Yerevan, Armenia 
'Florida International University, Miami, FL, USA 
8 Hampton University, Hampton, VA, USA 
9 University of the Witwatersrand, Johannesburg, South Africa 
10 University of Houston, Houston, TX, USA 
11 University of Virginia, Charlottesville, VA, USA 
12 Mississippi State University, Mississippi State, MS, USA 
13 James Madison University, Harrisonburg, VA, USA 

1Jt TRIUMF, Vancouver, British Columbia, Canada 
15 University of Regina, Regina, Saskatchewan, Canada 
16 Ohio University, Athens, OH, USA 
1 Tohoku University, Sendai, Japan 
18 Basel University, Basel, Switzerland 
19 Yerevan Physics Institute, Armenia 

20 Central Michigan University, Mount Pleasant, MI, USA 

21 College of William and Mary, Williamsburg, VA, USA 

22 Norfolk State University, Norfolk, VA, USA 
23 University of Manitoba, Winnipeg, Manitoba, Canada 
2i < Rensselear Polytechnic Institute, Troy, NY, USA 
(Dated: August 12, 2009) 

Cross sections for the p(e,e'7r + )n process on 1 H, 2 H, 12 C, 27 A1, 63 Cu and 197 Au targets were 
measured at the Thomas Jefferson National Accelerator Facility (Jefferson Lab) in order to extract 
the nuclear transparencies. Data were taken for four-momentum transfers ranging from Q 2 = l.l to 
4.8 GeV 2 for a fixed center of mass energy of W=2.1A GeV. The ratio of o\l and or was extracted 
from the measured cross sections for X H, 2 H, 12 C and 63 Cu targets at Q 2 = 2.15 and 4.0 GeV 2 
allowing for additional studies of the reaction mechanism. The experimental setup and the analysis 
of the data are described in detail including systematic studies needed to obtain the results. The 
results for the nuclear transparency and the differential cross sections as a function of the pion 
momentum at the different values of Q 2 are presented. Global features of the data are discussed 
and the data are compared with the results of model calculations for the p(e, e'n + )n reaction from 
nuclear targets. 

PACS numbers: 14.40.Aq,11.55. Jy,13.40.Gp,13.60.Le,25.30.Rw 



I. INTRODUCTION 

A fundamental challenge in nuclear physics is try- 
ing to understand the structure of hadrons in terms of 
their quark-gluon constituents, which are governed by 
the underlying theory of strong interaction called Quan- 
tum Chromodynamics (QCD). Measurements of exclu- 
sive processes such as pion clcctroproduction make it 
possible to extract meson form factors and study the 
quark-gluon distributions in the nucleoli. However, in or- 



der to develop a description of the atomic nucleus based 
on QCD, one also needs to understand how the prop- 
erties and interactions of hadrons change in the nuclear 
medium. Measurements of exclusive processes in the nu- 
clear medium sometimes help one to understand the in- 
teractions because QCD has definite predictions for ex- 
clusive processes in the medium. One such prediction is 
the phenomenon of color transparency (CT). 

In the study of CT, we can observe the interplay be- 
tween the creation of quark systems with small transverse 
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size, as we would expect from QCD, and the possible sup- 
pression of interactions of such point-like configurations 
(PLC) with the nuclear medium. A further complication 
is introduced by the formation time of the final hadron, 
which limits the lifetime of the PLC. To understand this 
phenomenon, it is thus important to study the dependen- 
cies on the four-momentum transfer squared of the vir- 
tual photon, Q 2 (size of the initial PLC, typically scales 
as ~ 1/Q), hadron momentum (formation length), and 
target mass A (path length through the medium). 

Due to the higher probability to form a PLC, a two- 
quark rather than a three-quark system is preferable, 
making pion electroproduction a natural choice for an 
initial study. However, pion electroproduction introduces 
the assumption of a quasi-free interaction of the photon 
with a virtual pion in the medium. The plausibility of 
this assumption can be addressed by comparing the ra- 
tio of the longitudinal to transverse cross sections from 
a nucleus with those obtained using a nucleon target, as 
differences in the behavior of this ratio would indicate 
effects of the nuclear medium resulting in the breakdown 
of the quasi-free assumption. 

Nuclear transparency is a natural observable in the 
quest to identify the onset of color transparency. Nuclear 
transparency is defined as the ratio of the cross section 
per nucleon for a process on a bound nucleon inside a 
nucleus to that from a free nucleon. Since the nuclear 
medium is not opaque to hadrons, color transparency will 
lead to an increase in nuclear transparency as a function 
of momentum transfer, hadron momentum and target 
mass. This simultaneous increase of nuclear transparency 
is distinct from those from other conventional processes 
such as rescattering. Moreover, it illustrates the need for 
careful experiment design that simultaneously measures 
the Q 2 and A dependence of nuclear transparency and at 
the same time tests alternative reaction mechanisms that 
can produce potential increases in nuclear transparency. 

With the availability of high-intensity, continuous elec- 
tron beams up to 6 GeV at Jefferson Lab it became pos- 
sible for the first time to determine simultaneously the A 
and Q 2 dependences of the differential pion cross sections 
for Q 2 in the 1-5 GeV 2 range. These data were acquired 
in 2004 for 1 H, 2 H, 12 C, 27 Al, 63 Cu and 197 Au targets. 
The A and Q 2 dependence of the extracted transparency 
were published in [![■ 

The purpose of this work is to describe the experi- 
ment and analysis and to present and discuss additional 
results. The article is organized as follows: the next sec- 
tions describe methods for extracting nuclear transparen- 
cies and discuss some earlier data. Section ITVl introduces 
the basic cross section formalism of pion electroproduc- 
tion. In section [V] the experiment is described including 
the experimental setup. The data analysis is described in 
section I VII The determination of the cross sections and 
the extraction of the nuclear transparencies and their as- 
sociated uncertainties are described in section I VIII and 
section IVIIII The results are presented in section IIXI 
The global features of the cross section ratio and the nu- 



clear transparencies are discussed and a comparison is 
made with results of recent calculations. 



II. METHODS OF DETERMINING NUCLEAR 
TRANSPARENCY FROM DATA 

In a simple picture of pion electroproduction from the 
nucleus, the electron exchanges a virtual photon with 
a proton which is moving due to its Fermi momentum. 
The struck proton turns into a neutron ejecting a pos- 
itively charged pion (quasi- free approximation). In the 
quasi-free picture, the ejected pion may interact with the 
residual nucleons and the fraction of pions which can es- 
cape from the nucleus is the nuclear transparency of pi- 
ons. In reality, the scattering process is more complicated 
and the deviations from this picture reveal much about 
the nucleus and their constituents. In the simple quasi- 
free picture, the ratio of the longitudinal to transverse 
cross section from a bound proton inside the nucleus is 
expected to be the same as that from a free proton. 

Assuming the dominance of the quasi-free process, one 
can extract the nuclear transparency of the pions by tak- 
ing the ratio of the acceptance corrected cross sections 
from the nuclear target to those from the proton. Nuclear 
transparency quantifies the interactions of the outgoing 
pion with the residual nucleons on its way out of the nu- 
clei, and is thus the key observable in searching for CT 
effects. Color Transparency was first proposed by Brod- 
sky and Mueller Q in 1982, and refers to the vanishing 
of the hadron-nucleon interaction for hadrons produced 
in exclusive processes. At sufficiently high momentum 
transfers, where the hadron is produced with small trans- 
verse size, b± f=a 1/Q 3], the fast-moving hadron also has 
a compact longitudinal size due to Lorentz contraction. 
Such a reduced size quark-gluon state is called a point- 
like configuration (PLC). The formation length, the dis- 
tance over which the PLC travels before the hadron 
reaches its final size can be written, 

j P ' ^lifetime /-. \ 

assuming a linear expansion in time, where f3 is the 
speed of the hadron in the lab frame and Tnf e ti me is the 
life time of the PLC in its rest frame. 

QCD predicts that the cross section for the interaction 
of a small qq dipolc is proportional to b 2 in the leading 
order approximation, where b is the transverse separation 
between the q and the q. 

If the qq dipolc remains small (pointlike) over the range 
of the nuclear system, the nucleus will be transparent 
to the produced hadron. This can be accomplished at 
large momentum transfers where the formation length is 
sufficiently large. 

At finite energies, however, the mechanism for the ex- 
pansion and contraction of the interacting small system 
is more complex. The quantum diffusion model [J] as- 
sumes that quarks separate in the transverse direction 
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at the speed of light and that their quark separation is 
proportional to y 7 ^, where z is the longitudinal distance 
from the production point to the position of the particle. 
The formation length in this model is determined from 
the average value of the dominant energy denominator. 

where Ph and Mh are the momentum and mass of the 
hadron, respectively, while M n is the mass of a typical 
intermediate state of the hadron. 

The precise value for AM 2 — (M 2 — M 2 ) -1 is a matter 
of great uncertainty 0, [|[ with estimates between 0.25 
and 1.4 GcV 2 . 

The onset of color transparency is expected at lower 
energies for production of mesons, for instance of pi- 
ons, than for production of baryons, because the quark- 
antiquark pair is more likely to form a small size ob- 
ject @. 

In pion clcctroproduction, the coherence length, de- 
fined as the distance over which the virtual photon fluc- 
tuates into a qq pair, is smaller than the nucleon radius 
at JLab energies, and for the present experiment it is 
essentially constant. This removes any coherence length 
dependence through t-channcl n — p exchange, in which 
the incident electron scatters from a virtual photon bring- 
ing it on-shell. Such a coherence length dependence can 
mimic a CT-like energy dependence. 

Color transparency has recently been put into the con- 
text of a QCD factorization theorem for longitudinally 
polarized photons in meson clcctroproduction. It pre- 
dicts that for sufficiently high values of Q 2 , at fixed xb 
and fixed momentum transfer to the nucleon, —t, the 
amplitude can be expressed in terms of a short-distance 
(hard) process, a distribution amplitude describing the 
final state meson, and Generalized Parton Distributions 
(GPDs) Q. The latter describe the long-distance (soft) 
physics. 

The factorization theorem has been proven at asymp- 
totically high Q 2 , but showing its validity at finite, high 
Q 2 requires stringent tests of all necessary conditions. 

The unambiguous observation of the onset of CT is 
a critical precondition for the validity of the factoriza- 
tion theorem. This is because where CT applies the out- 
going meson retains a small transverse size (inter-quark 
distance) while soft interactions like multiple gluon ex- 
change between the meson produced from the hard in- 
teraction and the baryon are highly suppressed. Factor- 
ization is thus rigorously not possible without CT. 

III. EARLY EXTRACTIONS OF NUCLEAR 
TRANSPARENCY 

The first experiments designed to search for color 
transparency used the 12 C(p, 2p) reaction and wereper- 
formed at Brookhaven National Laboratory (BNL) [H, [H, 
HoT |. The nuclear transparency was defined as the ratio 



of the elastic p — p cross section in the nucleus and hy- 
drogen. This ratio initially increased as a function of 
the beam energy and then decreased, with a peak near 
9 GeV. Though this behavior was not predicted by tra- 
ditional nuclear physics calculations, it is currently not 
attributed to CT. Instead it is associated with nuclear 
filtering [ll|, [T^] , which is related to the suppression of 
the long range components in the wave function, or to 
the threshold for charm resonance production [l3[ . 

The nuclear transparency was also measured using the 
A(e, e'p) reaction at Stanford Linear Accelerator Center 
(SLAC) [3 HH and at Jefferson Lab [H, H3- In this 
case, the nuclear transparency was defined as the ratio of 
the experimental yield and the one obtained from Plane 
Wave Impulse Approximation. The ratio was found to be 
energy independent from Q 2 ss 2 GeV 2 to 8.1 GcV 2 for 
deuterium, carbon, iron, and gold targets, thus indicating 
no significant CT effect. The absence of the CT effect in 
this channel has been interpreted as an indication that 
the proton formation length may only have been as large 
as internucleonic distances, rather than the size of the 
nucleus in these experiments [l8j . 

CT measurements using coherent and incoherent p° 
production have been performed at Fermilab [l9T ] and 
more recently at DESY |2lj. The nuclear transparency 
was defined as T = <ja/A(Jh- The results from Fermilab 
were parametrized with a function T = A™ -1 , and at a 
first glance, the observed positive slope of a as a func- 
tion of Q 2 appeared to contradict the flat Q 2 dependence 
predicted by the Glauber multiple scattering mechanism. 
However, the results have since been interpreted as a co- 
herence length effect [2(J, in which for kinematics with 
large coherence length, the virtual qq pair may undergo 
interactions with the nucleus before the hard interaction 
that puts it on-shell. The kinematics of the Fermilab 
experiment were not at a constant l c , and thus the vari- 
ation of a with Q 2 was attributed to a reduction in the 
initial state interactions rather than to a reduction in 
final state interactions. The more recent measurements 
from DESY, ranging between Q 2 =0.9-3 GeV 2 and at con- 
stant l c , showed a rise in nuclear transparency with Q 2 
consistent with theoretical calculations of CT [2~fl ] . 

The most convincing evidence for the existence of color 
transparency comes from an experiment performed at 
Fermilab [22|. There, the cross section of diffractive dis- 
sociation of 500 GeV/c pions into dijets was measured 
and parametrized with a=aoA a , where <jq is the tt — N 
cross section in free space. The free parameter, a, was 
fit to the data with the result a ~1.6. This result was in 
agreement with calculations assuming 100% color trans- 
parency and was very different from the normal n — N 
cross section, which has a dependence a—a^A 2 ^. 

At low energies, hints of CT effects have also 
been observed from pion photoproduction from helium, 
4 H e(7 , 7T~ p) , at Jefferson Lab j23[. CT can be measured 
in these reactions where Q 2 =0 by measuring the cross 
section vs. the four momentum transfer squared, t to 
the hadron system. The data showed 2a deviations from 
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where 



FIG. 1: Kinematics of the p(e,e' ,tv + )n reaction in the labo- 
ratory frame. 



the traditional Glauber calculations and the slope of the 
data vs. —t was in better agreement with calculations 
including CT. 



IV. ELEMENTARY PION 
ELECTROPRODUCTION 

A. Kinematics 

The kinematics of the pion electroproduction reaction 
studied here are shown in Figure [1] The incident elec- 
tron has a three momentum of k. The scattered electron 
has a three momentum k' and travels at a polar angle 
8 e in the laboratory frame with respect to the direction 
of the incident beam. The three momentum vectors of 
the incoming and outgoing electron define the scattering 
plane. The corresponding four-momenta are k = (E, k) 
and k! = (-E',k'). The virtual photon carries a four- 
momentum q = (u>, q), which is given by q = k — k' . The 
reaction plane is defined by the three-momentum q and 
the three-momentum vector of the produced pion P T and 
makes an angle 4>-n with respect to the scattering plane. 
The angle (in the lab system) between and q is 6 V . 

The pion electroproduction reaction can be described 
using three Lorentz invariants. In addition to Q 2 , we use 
the invariant mass of the virtual photon-nucleon system, 
W, which can be expressed as W = ^ / M 2 



2M p lj 



where M p is the proton mass, and t = {P^ — q) £ 
Mandelstam variable. 



is the 



B. Cross sections 

The pion electroproduction cross section from a sta- 
tionary proton in the one photon-exchange approxima- 
tion is 12411: 



d 5 <T 



d£l e >dE e idQ<T T 
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dfljr 
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a E e , K eq 1 
2tt 2 E e ~Q 2 1 - e 
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is the virtual photon flux, where a is the fine structure 



constant. The factor K e , 
equivalent photon energy, and 



{W 2 - M 2 )/(2M P ) is the 



tan^ — 
2 
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is the longitudinal polarization of the virtual photon. 
The two-fold differential cross section for a stationary 
proton target can be written as 



d 2 G 
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J- 
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where the solid angle of the pion, £1^ is determined in the 
lab frame, and J is the Jacobian for the transformation 
from to t, 4>. The two-fold cross section can be sepa- 
rated into four structure functions, which correspond to 
the polarization states of the virtual photon, a longitudi- 
nal (L), a trans verse (T), and two interference terms (LT 
and TT): 



2tt- 
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The interference terms vanish in parallel kinematics 
(6^ = 0) because of their dependence on V (25[. 

The four structure functions can be separated if mea- 
surements at different values of e and </>„■ are performed 
(L/T separation), where W,Q 2 and t are kept constant. 
The photon polarization e can be varied by changing the 
electron energy and scattering angle. 
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1.10 


2.26 


0.050 


4.021 


27.76 


1.190 


10.61 


2.793 


0.23 


0.50 


2.15 


2.21 


0.158 


5.012 


28.85 


1.730 


13.44 


3.187 


0.41 


0.56 


3.00 


2.14 


0.289 


5.012 


37.77 


1.430 


12.74 


3.418 


0.56 


0.45 


3.91 


2.26 


0.413 


5.767 


40.38 


1.423 


11.53 


4.077 


0.70 


0.39 


4.69 


2.25 


0.527 


5.767 


52.67 


1.034 


10.63 


4.412 


0.79 


0.26 


2.16" 


2.21 


0.164 


4.021 


50.76 


0.730 


10.60 


3.187 


0.42 


0.27 


4.01 6 


2.14 


0.441 


5.012 


55.88 


0.910 


10.55 


3.857 


0.71 


0.25 


2.16 c 


1.74 


0.374 


4.021 


32.32 


1.730 


19.99 


2.074 


0.65 


0.63 



"low e setting used for L-T separations 
''low e setting used for L-T separations 
c test point to check the dependence on W and 

TABLE I: The central kinematic settings for the pion trans- 
parency experiment, q is the angle between the virtual photon 
3-momentum q and the beam direction in the lab frame, 
is the magnitude of the three momentum of the virtual struck 
pion in the quasi-free knockout picture, and x — Q 2 / (2Mv) . 
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For nuclei, where there is a new degree of freedom due 
to the Fermi momentum of the struck nucleon, the dif- 
ferential pion clcctroproduction cross section is given by, 



r„ 



d 3 a 



dVl e >dE e > dflirdPir " dQndPrr ' 



(8) 



can be 



The three-fold differential cross section, ,„ ° 

transverse and interference 



separated into longitudinal, 
terms just as in Eqns. 16171 



Quantity 


HMS 


SOS 


Max. Central Mom. 


7.3 GeV/c a 


1.74 GeV/c 


Optical Length 


26.0 m 


7.4 m 


Momentum Acceptance 


±10% 


±20% 


Solid Angle 6 


6.7 msr 


7.5 msr 


In-Plane Ang. Acc. a 


±27.5 mrad 


±57.5 mrad 


Out-of-Plane Ang. Acc.° 


±70 mrad 


±37.5 mrad 



"To date, the HMS has been operated to 6.0 GeV/c 
'The solid angle and angular acceptances are given for the large 
collimators in both the HMS and SOS spectrometers. 



EXPERIMENT AND SETUP 



TABLE II: Selected properties for the HMS and SOS. 



The E01-107 experiment (pionCT) was carried out in 
Hall C at Jefferson Lab [28[ in 2004. A schematic view of 
the experimental setup is shown in Fig. [2] and the kine- 
matic settings of the measurements are shown in Table [I] 
The continuous wave (100% duty factor) electron beam 
from the CEBAF accelerator has a characteristic 2 ns 
micro-structure that arises from the 1.5GHz rf structure 
of the accelerator and the 3 hall beam splitting scheme. 
This 2 ns structure is used in the analysis to identify 
(e, e'ir + ) coincident events. The electron beam with en- 
ergies of up to 5.8 GeV was incident on liquid hydro- 
gen and deuterium, and solid foil targets of 12 C, 27 A1, 
63 Cu and 197 Au. For the cryotargets, a 4.0 cm diameter 
cylindrical cell with an axis perpendicular to the beam 
direction was used. The cell walls were made from an alu- 
minum alloy with a thickness of 0.01 cm. The integrated 
beam current (35-80 ^A) was measured using a pair of 
microwave cavities with a DC current transformer to an 
accuracy of 0.4%. To reduce target density fluctuations 
in the cryotargets arising from beam heating, the beam 
was rastered to a 2x2 mm 2 profile. No target density 
fluctuations were measured to an accuracy of 0.6%. 

The scattered electrons were detected in the short 
orbit spectrometer (SOS). The SOS consists of three 
non-superconducting magnets, one quadrupole followed 
by two dipoles, which share a common yoke. The 
quadrupole focuses in the non-dispersive direction, the 
first dipole bends particles with the central momentum 
upwards by 33°, and the second one bends them down- 
wards by 15°. The SOS spectrometer design was opti- 
mized for detection of unstable and short-lived particles 
and thus has a flight path of only 10 m. The electropro- 
duced pions were detected in the high momentum spec- 
trometer (HMS). The HMS consists of three supercon- 
ducting quadrupole magnets followed by a momentum 
analyzing vertical-bend superconducting dipole used in 
a point-to-point tune for the central ray. The detector 
package is mounted near the focal plane of the optical 
system, which is located inside a concrete shield house. 
The HMS has a 26 m pathlcngth and a maximum central 
momentum of 7.3 GeV/c A detailed description of the 
spectrometers and the spectrometer hardware is given 
in [2^, [3(| • Selected properties of the two spectrometers 
are listed in Table [TTJ 



Both spectrometers are equipped with multiwire drift 
chambers for particle tracking, and segmented scintilla- 
tor hodoscope arrays for time-of-flight measurements and 
triggering. 

The HMS has a gas Cerenkov and a lead glass calorime- 
ter in the detector stack for K + /ir + separation. For this 
experiment, the Cerenkov was filled with C4F10 gas at 
97 kPa. The corresponding index of refraction is 1.00137 
resulting in momentum thresholds of 2.65 GeV/c for ir + 
and 9.4 GeV/c for K + . An aerogel Cerenkov detector [3l[ 
was also used in the detector stack. It was used to iden- 
tify 7r + for central momentum settings below 3.1 GeV/c. 
The aerogel had an index of refraction of 1.015 giving 
thresholds of 0.8 GeV/c for pions and 2.85 GeV/c for 
kaons. 

The SOS has a combination of a gas Cerenkov and a 
lead-glass calorimeter for e~ /ir~ separation. The nomi- 
nal SOS Cerenkov detector is filled with CCI2F2 gas at 
101 kPa with an index of refraction of 1.00108. Dur- 
ing part of the experiment (E e =5.767 GeV) the nom- 
inal SOS Cerenkov detector was replaced by one filled 
with C4F10 gas at 143 kPa with an index of refraction of 
1.0019. The corresponding thresholds for electrons were 
below 10 MeV/c and those for pions were 3 GeV/c (nom- 
inal) or 2.27 GeV/c (new). 



VI. DATA ANALYSIS 

The raw data were processed and combined with addi- 
tional information from the experiment like momentum 
and angle settings of the spectrometers, detector posi- 
tions, and beam energy to give the particle trajectories, 
momenta, velocities, energy deposition, and to perform 
particle identification. 



A. Event Reconstruction 

The spectrometer quantities x, y, x' and y' ', were de- 
duced from reconstruction of the wire chamber data. 
These quantities are the vertical and horizontal positions 
of the track at the midpoint between the wire chambers, 
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FIG. 2: Schematic view of the experimental setup for E01-107. 



and the gradients of that track with respect to the spec- 
trometer central ray. The target quantities ytar, x' tar , 
y' tar and 5 were determined from the spectrometer quan- 
tities by suitable transformation functions. These quan- 
tities are the horizontal position of the event, the hor- 
izontal and vertical gradients of the track with respect 
to the spectrometer central ray and the momentum of 
the particle given as a percentage above the central mo- 
mentum setting of the spectrometer, respectively. A spe- 
cial data set, using a series of foils targets placed at well 
determined positions in the target region was taken to 
calibrate the transformation matrix. 



B. Particle Identification and event selection 



rejection ratio for the HMS Aerogel Cerenkov detector 
was approximately 5:1 at P cen trai = 2.1 and 2.8 GeV/c. 
A cut on the gas Cerenkov detector varying between 

Nphotoelectrons >0.7 and ~N p hotoelectrons >2.0 was USed 

when Prms > 3.2 GeV/c. The corresponding cut ef- 
ficiencies were determined using tight cuts on the coinci- 
dence time and the aerogel Cerenkov detector to remove 
protons, and were found to be > 98.2±0.5% The cut ef- 
ficiency was parametrized as a function of the HMS frac- 
tional momentum, Shms to take into account that pions 
at negative Shms are closer to the momentum threshold 
of the detector compared to those at positive Shms- The 
resulting rejection ratio was 50:1 at P ce „ tra ;=3.2 GeV/c 
and 300:1 at P ce ntrai = 4.4 GeV/c with an uncertainty 
of 0.2%. 



Electrons were selected with a SOS gas Cerenkov cut 

Of Nphotoelectrons > 1.0 for the nominal Or Nphotoelectrons 

> 5.0 for the new detector. The efficiency of the cut 
was determined, for each of the SOS central momentum 
settings, using a sample of electrons that was identified 
using the lead-glass calorimeter. The resulting efficiency 
was always higher than 99.2% with an uncertainty of 
0.2%. The pion rejection ratio was 100:1 for the nom- 
inal and 300:1 for the new detector. In the HMS, the 
aerogel and gas Cerenkov detectors were used to select 
7r + . The aerogel was used for additional particle identifi- 
cation when the central momentum setting of the HMS, 
Phms, was less than 3.2 GeV/c. The lower limit of 
the aerogel efficiency was determined using tight cuts 
on the coincidence time, and was found to be >98.8± 
0.5% for a threshold cut of Nphotoeiectrons >0.7. The 



C. Efficiencies 

The raw yield was normalized by the beam charge, 
and the data were corrected for detector efficiencies, com- 
puter dead time, correction to target thickness and nu- 
clear absorption in the detection materials. 

1. Tracking efficiency 

The tracking efficiency is the probability of finding a 
track from experimental signals from the wire chambers 
when a charged particle passes through them. It de- 
pends on the efficiency of the wire chambers and also 
on the algorithm used to construct a track for a given 
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event. The tracking efficiency was >92% for the HMS 
and >96% for the SOS, and was calculated run-by-run 
for this experiment. The uncertainties to the tracking 
efficiencies are estimated from the stability of the track- 
ing efficiency as a function of the total rate in the spec- 
trometer, and were found to be 1.0 % for the HMS and 
0.5% for the SOS. The rate dependence of the tracking 
efficiency, which is related to the increased probability to 
have multiple tracks at high rates was taken into account 
with a tracking efficiency calculation including multiple 
track events [MlH- 

2. Trigger Efficiency 

The HMS and the SOS spectrometers each contain four 
layers of scintillators and a requirement for a single-arm 
trigger is that 3 out of 4 layers must have a hit. The 
3-out-of-4 efficiency was found to be above 99.5% for all 
runs in both spectrometers with an uncertainty of 0.5% 
for the HMS and 0.5% for the SOS. 



3. Computer and Electronics dead times 

The computer dead time, which are due to the fi- 
nite time required by the computers to process an event, 
can be directly calculated from the number of generated 
pretriggers and the number of accepted triggers. The 
computer dead time during the pionCT experiment was 
about 25%. The electronics dead time, due to the finite 
time required by the electronic modules was estimated 
from copies of the original pretrigger signal at varying 
gate widths. The correction due to electronics dead time 
was less than 1%. The total uncertainty (0.2-0.5%) due 
to the dead time corrections was dominated by the com- 
puter dead time. 

4. Coincidence Blocking 

Coincidence events can be blocked when a non- 
coincident event arrives in one of the spectrometers just 
before the coincident event. The "coincidence blocking" 
events are lost from the data due to the coincidence time 
cuts in the analysis. The coincidence blocking correction 
was estimated from the rate dependence of the number 
of blocked events, and was found to be less than 0.7% 
with an uncertainty of about 0.2%. 

5. Pion Absorption 

Pions may interact through the strong nuclear force 
with the nucleus, for instance in the target material, the 
window of the scattering chamber, or the windows of the 
spectrometer. These events are lost before reaching the 
detectors in the HMS detector hut. We account for these 



events by correcting for pion absorption in the HMS. The 
transmission for pions in the HMS is about 95% [2(| and 
depends weakly on the pion momentum between 2.1 and 
4.4 GeV/c. In addition some pions interact with the 
scintillator material producing events that reconstruct 
at lower pion velocity {(3 n ) or even 0^ = and lower 
coincidence time. Using a 2-dimensional cut on (3 V and 
coincidence time, such events were included in our pion 
yield. The efficiency of such a cut was studied using 
(e, e'ir~) [32| and H(e,e'p) [H, HU data. The difference 
between the 2-dimcnsional cut and a simple cut on pion 
velocity that does not include the pion with lower velocity 
and the (3^ =0 events was found to be within the uncer- 
tainty associated with the HMS detector efficiency and 
the uncertainty of the pion absorption correction. The 
uncertainty in the pion absorption correction was esti- 
mated from the difference between the calculated pion 
transmission and the measured proton transmission, and 
was found to be 2% in the absolute cross section determi- 
nation [2(| Hi[ . The difference in target thickness leads 
to an additional A dependent uncertainty which is esti- 
mated to be 0.5%. 



6. Backgrounds 

The coincidence time was calculated from the time dif- 
ference between the SOS and HMS triggers, and was 
used to help identify e~ /tt + coincidences from e~ /p co- 
incidences. Corrections to the coincidence time include 
the path length of the tracks through the magnetic el- 
ements of the spectrometer, difference in signal propa- 
gation times in cables, pulse height corrections for the 
signals from the scintillators, and subtraction of the time 
required for the light to travel in the scintillators from 
the event position to the PMT. Random coincidences, 
resulting from an electron and pion from different beam 
buckets, have a coincidence time that is offset from the 
in-time by multiples of 2 ns. The in-time e~ /ir + events 
were selected with a cut on the coincidence time around 
the central peak. The random background was estimated 
by averaging over three bunches, that could be readily 
identified by eye, to the right and three bunches to the 
left of the in-time peak, and then subtracting it from the 
in-time yield. 

Background from the cryotarget cell walls, less than a 
few percent, was measured and subtracted using an alu- 
minum target of approximately seven times the thickness 
of the target cell walls. The contribution of the cell walls 
was small (less than 5%), and due to the high statistical 
accuracy of the dummy target data, the contribution of 
the subtraction was < 0.1% to the total uncertainty. 

7. Missing mass 

Once true e~ /ir + coincidences are identified, the miss- 
ing mass of the recoiling nucleoli system was recon- 
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structed from the measured quantities. In the present 
analysis, a cut in the missing mass was used to ensure 
that no additional pions were produced in the case of hy- 
drogen, and to minimize the contribution of multi-pion 
events to less than 5% with an uncertainty of < 0.4% (see 
section IVII A 41 for a more detailed discussion). 

VII. DETERMINATION OF THE NUCLEAR 
TRANSPARENCY 

To extract cross section information from the data, the 
measured yields were compared to the results of a Monte 
Carlo simulation for the actual experimental setup (see 
section IVII A|) , which included a realistic model of the 
pion electroproduction cross section. When the model 
describes the dependence on W, Q 2 , t, and 9 V (P-k) of 
the four structure function in equation [7] well, the cross 
section for any value of W and Q 2 in the acceptance can 
be determined as 

(da\ = Ydatafda\ 
\dtj data Y mc \dt) mc 

where Y is the yield over W and Q 2 . The term "data" 
refers to the measured experimental yield, and the term 
"mc" refers to the simulated events and yield. 

To extract the nuclear transparency, which is defined 
as the ratio of cross sections extracted from data and 
from a model of pion electroproduction from the nucleus 
without 7r — N final state interactions, simulations were 
performed as described in the next section. 

A. Monte Carlo Simulation 

The standard Monte Carlo simulation code for Hall 
C, SIMC, was used to simulate the experiment. Events 
were generated in a phase space marginally larger than 
the acceptance of the spectrometers. After events were 
generated at the vertex, they were transported through 
the spectrometer optics using COSY matrix elements de- 
termined from a COSY INFINITY [35] model of the 
spectrometers. Each event was weighted by the relevant 
model cross section (see section IVII A 2|) . A comparison 
between experimental and simulated distributions of re- 
constructed quantities for a hydrogen target are shown in 
Fig. [31 If the detector setup and the spectrometer accep- 
tances (including coincidence acceptance) is realistically 
modeled in a simulation, the boundaries of the distri- 
butions should match. Differences in magnitude can be 
attributed to differences between the actual cross section 
and the one used in the model. 

To describe electroproduction from nuclear targets the 
quasi-free approximation was used, since the energy of 
the incoming electron is large compared to the energy 
associated with the binding of the nucleons. Properties 
of the nucleons inside the nucleus are assumed to be de- 
scribed by an independent particle shell model, where 
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FIG. 3: (color online) Experimental (crosses) and Monte 
Carlo (lines) distributions for the hydrogen target at Q 2 =3.91 
GeV 2 . 



each nucleon interacts with a mean field exerted by the 
other nucleons. The probability of finding a nucleon with 
momentum p m , and separation energy E m , in the nu- 
cleus is given by a spectral function, S(E m ,p m ). Exist- 
ing s pec tral functions were used for deuterium, carbon, 
gold [l4 [lU , and aluminum [36| . The copper spectral 
function was constructed from the iron spectral function 
described in Ref . by increasing the number of protons 
in the outermost 1/ shell from 6 to 9 and changing the 
central binding energy of this shell using the separation 
energy for copper. The spectral functions did not include 
any corrections to account for shifts in their strength to 
large missing momentum caused by nucleon-nucleon cor- 
relations. 

Although the momentum of the proton is given by p m , 
E m is not constrained by any of the assumptions in the 
quasi-free approximation. In the present analysis, the 
off-shcllness was described by, 

M A = M p + M\_ x 

E P = AI A - ^(MI-i) 2 + |Pm| 2 (10) 

where Ma is the nuclear missing mass, M p is the mass of 
the nucleon, M* A _ X is the mass of the spectator nucleons, 
and E p is the energy of the struck proton. A comparison 
between experimental and simulated distributions of re- 
constructed quantities for nuclear targets is shown in in 
Figures H and 

The radiative corrections used in this analysis are 
based on the formalism of Ref. [38[ , which was explicitly 
modified and improved for (e, e'p) coincidence reactions 
in Ref. [3{|. It includes both external and internal radi- 
ation. The formalism was further modified for pion elec- 
troproduction, where the target particle is assumed to be 
a stationary proton, while the final pion is assumed to be 



9 



off-shell. The assumption of off-shellness was tested in 
the simulation and was found to change the Monte Carlo 
yield by at most 0.5%. The largest source of uncertainty 
in the simulation of radiative processes comes from radia- 
tion due to the pion as the electron radiation is relatively 
well known. The simulated yield changed by 2-4% when 
radiation from the pion was turned off depending on A 
and Q 2 . For the present analysis, we assume a global 
uncertainty due to radiative correction of 2%. The un- 
corrected uncertainty in the radiative corrections was 
estimated from the target dependence of the simulated 
yield when the pion radiation was turned off, and was 
found to be 1% at low Q 2 and 2% at high Q 2 . 

SIMC incorporates the effects of the pion decay, multi- 
ple scattering and energy loss. For the present analysis, 
2.5% (1.4%) of pion decays occur inside a magnet at the 
lowest (highest) Q 2 setting. The random uncertainty is 
about 0.5% and mainly accounts for muons coming from 
pions normally outside the acceptance. The difference 
between targets is smaller than 0.1%. The Monte Carlo 
also takes into account pions which punched through the 
spectrometer collimators. The simulation is based on the 
calculation of the pion transmission through materials 
described in Ref. [341 ] . 

For the present analysis, SIMC was modified to take 
into account the effect of Pauli blocking. This correction 
which is determined using the distribution function cal- 
culated by Fantoni and Pandharipande [4(| is applied as 
a weight to each event. 

The final state interactions between the knocked-out 
neutron and the residual nucleons (n-N FSI) can effect 
the quasi-free cross section and can shift strength in the 
missing mass spectrum near the single pion production 
threshold. It is likely to have the strongest effect when 
the relative momentum between the recoiling neutron 
and the spectator nucleons is small. An earlier exper- 
iment on light nuclei at JLab has shown that the effect 
of n-N FSI on the quasi-free cross-section reduces with 
increasing Q 2 (4l| . Since this experiment was conducted 
at relatively high Q 2 the n-N FSI was not explicitly ac- 
counted for in the simulation. One can account for a 
fraction of the shift in strength in the missing mass spec- 
trum caused by n-N FSI, by setting the missing energy 
E m = 0, when calculating the energy of the proton. How- 
ever, the lack of full accounting of the n-N FSI is most 
likely to be the cause of the discrepancy in the shape of 
the missing mass spectrum at the lowest Q 2 = 1.1 GeV 2 
setting, seen in Fig. [5] and 

Coulomb corrections to the incoming and scattered 
electron are applied according to the Effective Momen- 
tum Approximation (EMA) approach [42j. This ap- 
proach includes an improvement over earlier versions of 
the EMA, using an average potential to account for the 
focusing of the incoming electron wave function [42| . No 
Coulomb corrections were applied to hydrogen and deu- 
terium, since this effect is already included in the elemen- 
tary pion cross section. The corrections for the copper 
target ranges between 0.2-2%, while the yield correction 
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for the gold target ranges between 0.9-4.4%. We assume 
that 25% of this correction contributes to the uncorre- 
cted uncertainty 



1. Detector acceptances in SIMC 

The acceptances for both spectrometers used in this 
analysis are within the previously determined safe bound- 
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FIG. 6: (color online) Experimental (crosses) and Monte 
Carlo (lines) distributions for the copper target at Q 2 =3.91 
GeV 2 . The vertical line in the bottom right panel shows 
the position of the two-pion production missing mass cut 
(see \ VII A J\ for discussion). 



2. The model cross section 

The Monte Carlo was used to extract the bin-centered 
experimental cross section by iterating the model cross 
section until Monte Carlo distributions matched the data. 
The starting model was based on the previous pion elec- 
troproduction data from Hall C (34[. The model cross 
section was taken as the product of a global function 
describing the VF-dependence times (a sum of) Q 2 and 
t and dependent functions for the different structure 
functions. A correction function, which was assumed to 
factorize, 



C H {W,Q 2 ,t, 



0(W)K(Q 2 )T(t)F( ( / )pq ), (11) 



was determined by iterating the model and comparing it 
to the hydrogen elementary cross section. These correc- 
tion functions were assumed to be second order polyno- 
mials, with the exception of F((f) pq ), which was assumed 
to be of third order. 
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aries [34|. The phase space was studied by comparing 
distributions for reconstructed quantities like Q 2 , t, and 
missing mass as shown in Figure [3J 

The acceptance uncertainties were estimated using the 
over constrained H(e, e'p) reaction. These uncertain- 
ties arise from the uncertainty in the knowledge of the 
momentum and angle setting of the spectrometers and 
the beam energy. The point-to-point yield variation was 
found to be 1% for different Q 2 values. 



3. Iteration procedure 

For the hydrogen target the bin centered experimental 
cross section is given by 
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(12) 

For nuclear targets, the bin centered experimental cross 
section is given by 



dQ. e ' dE e i dVl^dPTr 



Y, 



data 



d 6 



&A 



Ymc KdSle'dEe/dCl^dPTr 



Xq 

(13) 

where the subscript, xq, indicates that the cross section 
is evaluated at a particular point (Wo, Qq, Oq) inside the 
acceptance. The model cross sections for hydrogen and 
for the nuclear targets were determined using a point 
Monte Carlo simulation, which was performed using scat- 
tered electron kinematics and pion angles generated ran- 
domly within a very narrow phase space volume that 
corresponds to Wo, Q\ and #o- For nuclear targets, 
is generated over the whole phase space from which a 
narrow range was selected in the analysis. 

The extracted cross sections depend on the initial cross 
section model, and thus there are systematic uncertain- 
ties associated with it. This uncertainty is obtained 
by extracting the cross section using a different starting 
model and was found to be 1.1%. However, this uncer- 
tainty will not contribute to nuclear transparency, since 
nuclear transparency involves the ratio of two Monte 
Carlo yields. 
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FIG. 8: (color online) Nuclear missing mass distributions (in 
GeV 2 ) for 12 C(e,e'-K + ) The data (red crosses) are compared 
to the simulation (blue line) which is a sum of single-pion and 
multiple-pion simulations. The shaded area (green) shows the 
contribution from the multi-pion simulation. The full simula- 
tion is normalized to the data. The dashed vertical lines repre- 
sents the threshold for double-pion production (11.34 GeV 2 ). 
The solid lines represent the position of the cut used in this 
analysis. 



FIG. 9: (color online) Nuclear missing mass distributions (in 
GeV 2 )for 63 Cu(e,e'n + ). The data (red crosses) are compared 
to the simulation (blue line) which is a sum of single-pion and 
multiple-pion simulations. The shaded area (green) shows the 
contribution from the multi-pion simulation. The full simula- 
tion is normalized to the data. The dashed vertical lines repre- 
sents the threshold for double-pion production (59.33 GeV 2 ). 
The solid lines represent the position of the cut used in this 
analysis. 



4- Multi pion production in nuclear targets 

The quasi-free pion production model is limited to sin- 
gle pion production. The production of more than one 
pion in a single event (multiple pion production) was 
suppressed for hydrogen targets during the pionCT ex- 
periment due to the relatively high Q 2 > 1 GeV 2 and 
W >2.1 GeV being above the resonance region. This 
suggests that the mechanism for multiple-pion produc- 
tion involves the outgoing pion producing one or more 
pions from a nucleon in a second process that was inco- 
herent from the production of the first pion. Multiple- 
pion events can only be produced above a missing mass 
threshold that is larger than the missing mass threshold 
for single-pion production, i.e., M x = Ma-i + M n for 
a nucleus of mass A. Indeed, in the analysis, at first a 
missing mass cut exactly corresponding to this thresh- 
old was used to suppress multiple-pion events. However, 
this cut resulted in an unacceptable loss of statistics at 
the highest Q 2 settings, and thus an alternative cut on 
the nuclear missing mass above the multi-pion threshold 
were used. 

In order to describe events above the two pion thresh- 
old, a multiple-pion production simulation was developed 
for the present nuclear target analysis. The mechanism 
for multiple-pion production was assumed to be quasi- 
free single-pion production from a nucleon followed by 



a secondary process that was incoherent from the first, 
where the pion produces one or more pions from a differ- 
ent nucleon. The cross section for the second process was 
assumed to be uniform over the acceptance of the HMS 
spectrometer. 

The effect of multi-pion production can be seen in 
Fig. [5] and [5] for the carbon and copper targets. The 
agreement between the missing mass distributions ob- 
tained from data and simulation improves with increasing 
Q 2 . The discrepancy seen at Q 2 = 1.1 GeV 2 is attributed 
to the reaction mechanisms missing from the simulation 
such as final state interactions between the knocked-out 
neutron and the residual nucleons (nN-FSI) and short 
range correlations. The effect of these reaction mecha- 
nisms reduces with increasing Q 2 . 

These results show that it is safe to increase the double- 
pion missing mass cut above the threshold with minimal 
contamination. The double-pion missing mass cut was 
placed at the position where the systematic uncertainty 
from the contribution of multiple-pion events was less 
than 5%. With these cuts, the total uncertainty due 
to multi-pion contamination is < 0.4%. We also noted 
an interesting smooth A dependence in the ratio of the 
multiple-pion to single-pion yields. 
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FIG. 10: (color online) The plot on the upper left panel is 
the phase space comparison between the high (black) and low 
(blue) e for —t vs Q 2 . The plot on the upper right panel is the 
phase space comparison for W vs Q 2 . The plot on the lower 
left panel shows the comparison for —t vs W . The lower right 
panel shows the azimuthal coverage for these two data sets. 
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FIG. 11: (color online) Representative plot of the experimen- 

,2 

tal cross sections, dt ^ as a function of the azimuthal angle 

<j> at Q 2 =2.15 GeV 2 for high and low e. The curves shown 
represent the fit of the measured values of the cross section to 
equation [?[ The —t in this plot corresponds to the common 
region between the high and low e data shown in Fig. \10\ Only 
the statistical uncertainties are shown. 



5. Test of the quasi-free assumption 

The average cross sections were extracted by integrat- 
ing over the whole acceptance (W, Q 2 and t). This aver- 
aging reduces the systematic uncertainties related to the 
cross section model in the Monte Carlo by smearing the 
exact kinematic information of the extracted cross sec- 
tion. We bin the data in <f> by integrating over all other 
kinematic variables at each e setting. The azimuthal an- 
gular coverage for the hydrogen target at Q 2 = 2.15 GeV 2 
can be found in Fig. [TU1 Due to the correlation in the 
kinematics, the central W, Q 2 and t values for each 4> 
bin are different from those obtained by integrating over 
the entire 4> region. Thus, a Monte Carlo simulation with 
model cross sections is used to account for the correction 
between the cross section evaluates at the center of each 
</> bin and the one for the entire 4> region. The four struc- 
ture functions in equation [7] are extracted by fitting the 
data with respect to (f> for both high and low e settings 
simultaneously. A representative fit for the hydrogen tar- 
get at Q 2 = 2.15 GeV 2 can be found in Fig. [TH In the 
analysis, an additional acceptance cut is used to ensure 
that the kinematic region given by IF, Q 2 and t is the 
same at high and low e. Such a phase space comparison 
is shown in Fig. [TO] 

The same fitting procedure described above could be 
used to obtain the Roscnbluth-separated pion electropro- 
duction cross sections. However, where this separation is 
relatively straightforward for hydrogen targets, a similar 
separation for nuclear targets relies on the assumed quasi- 
free reaction mechanism. This is because beyond the W, 
Q 2 and t dependence, the elementary off-shell pion elec- 
troproduction cross section has a P^ dependence in the 
nuclear medium. Hence the input model for elementary 
pion electroproduction needs to starting value for the P w 
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FIG. 12: (color online) The ratio of longitudinal to transverse 
cross sections for pion electroproduction from H, 2 H, 12 C, 
and 63 Cu targets, at fixed Q 2 = 2.15 (left) and 3.91 (right) 
GeV 2 . The inner error bars represent the statistical uncer- 
tainties, while the outer error bars are the sum in quadrature 
of the statistical and systematic uncertainties. The curves 
represent constant-value fits to all nuclear data at a fixed Q 2 . 
The probability of these constant-value fits assuming Gaussian 
statistics is 69 and 70%, respectively. 



dependence. This quasi-free P^-dependence is taken into 
account in the Monte Carlo simulations for the nuclear 
targets, and the iterative procedure described earlier is 
followed. This implies that the extracted nuclear cross 
sections represent the averaged values integrated over a 
wide kinematic acceptance, rather than the bin centered 
value. These averaged cross sections are used to obtain 
the longitudinal to transverse ratio between nuclear and 
hydrogen targets. The ratio of longitudinal to trans- 
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Item 


Uncorrelated uncertainty 


A dependent 


Correlated uncertainty 




(%) 


(%) 


(%) 


HMS Cerenkov 


0.2 




0.3-0.5 


SOS Cerenkov 


0.2 




0.3-0.5 


Charge 


0.4-0.9 




0.4 


Coincidence blocking 


0.2 






HMS Trigger 


0.5 






Dead time 


0.2-0.5 






HMS Tracking 


1.0 




1.0 


SOS Trigger 


0.5 






SOS Tracking 


0.5 




0.5 


Pion decay 


0.1 




1.0 


Coulomb corrections 


<1.0 






Radiative corrections 


1.0-2.0 




2.0 


Collimator 


0.5 




1.0 


Acceptance 


1.0 




2.0 


Iteration procedure 


1.1 






Multi-pion contamination 


<0.4 






Target thickness 




0.5 - 1.0 




Pion absorption 




0.5 


2.0 


Total 


2.4-3.4 


0.7-1.1 


3.9-4.0 


model dependence 


3.5-7.6 







TABLE III: The systematic uncertainties in extracting nuclear transparency. The uncorrelated uncertainty contributes directly 
to the point-to-point uncertainty in the nuclear transparency, the A dependent uncertainty is independent of Q 2 , while the 
correlated uncertainty is independent of both the target nucleus and Q 2 . 



Item 


Uncertainty in 
Differential cross section 
(%) 


Uncertainty in 
L-T cross section ratio 
(%) 


HMS Cerenkov 


0.4-0.5 




SOS Cerenkov 


0.4-0.5 




Charge 


0.6-1.0 




Coincidence blocking 


0.2 




HMS Trigger 


0.3-0.5 




Dead time 


0.2-0.5 




HMS Tracking 


1.1-1.4 




SOS Trigger 


0.3 




SOS Tracking 


0.5 




Pion decay 


1.0 




Coulomb corrections 


<1.0 




Radiative corrections 


2.2-2.8 




Collimator 


1.1 




Acceptance 


2.2 




Iteration procedure 


1.3 - 1.5 


13.0-18.0 


Multi-pion contamination 


<0.4 




Target thickness 


0.5-1.0 




Pion absorption 


2.1 




Kinematics 


1.5-2.0 


3.0-16.6 


momentum coverage 




<12.0 


Total 


4.8-5.7 





TABLE IV: The systematic uncertainties in extracting cross sections and in extracting the ratio of longitudinal to transverse 
cross section 
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verse cross sections at fixed Q 2 = 2.15 and 3.91 GeV 2 
for the various targets used in this experiment are shown 
in Fig. 1121 We find no difference, within the experimental 
uncertainties, between the ratio of longitudinal to trans- 
verse cross sections for nuclear and hydrogen targets, this 
can be viewed as a confirmation of the quasi-free reaction 
reaction mechanism. 



VIII. UNCERTAINTY ESTIMATES 

Table Mil lists the systematic uncertainties in ex- 
tracting the nuclear transparency. Several sources of 
these uncertainties have already been discussed in sec- 
tions EnXTEnXl 

The uncertainty in the acceptance is based on exten- 
sive single-arm elastic and deep- inelastic measurements 
from [H, and 1 H(e,e'p) data, including sieve-slit 
data on a carbon target, taken to check the optical ma- 
trix elements. The influence of the uncertainties in the 
offsets in the kinematical variables such as beam energy, 
momenta and angles, were determined by changing their 
values by their uncertainty and evaluating the resultant 
changes in the cross sections. 

The uncertainty in the target thickness for the solid 
targets is dominated by the knowledge of their purities 
and the thicknesses. For the 4 cm liquid targets, the un- 
certainty is dominated by the target boiling correction 
and the beam path length since the beam did not pass 
through the axis of the liquid target cylinder. The un- 
certainty in the total charge includes a 0.4% uncertainty 
due to calibration of the beam current monitors. 

The largest correlated systematic uncertainties are ra- 
diative corrections, and acceptance, resulting in a total 
correlated uncertainty of 3.9-4.0%. The uncorrelated sys- 
tematic uncertainty is dominated by contributions from 
the acceptance, radiative corrections and the iteration 
procedure, resulting in a total uncorrelated uncertainty 
of 2.4 - 3.4%. The largest contribution to the " A depen- 
dent" uncertainty is the target thickness, resulting in a 
total uncertainty of 0.7 to 1.1%. 

The largest source of uncertainty in the radiative cor- 
rection procedure comes from pion radiation as the elec- 
tron radiation is relatively well known. The Monte Carlo 
equivalent yields changed by 2-4% when pion radiation 
was turned off (2% at low Q 2 and 4% for the heavy tar- 
gets at high Q 2 ). From this the normalization uncer- 
tainty was taken to be 2%. The point-to-point uncer- 
tainty in the radiative corrections was estimated from the 
target dependence of the Monte Carlo equivalent yield 
when the pion radiation was turned off. This was 1% at 
the low Q 2 and 2% at high Q 2 . 

In addition to these uncertainties, a Q 2 dependent 
model uncertainty was also determined to be 3.5% -7.6%. 
This uncertainty is the quadrature sum of the change 
in Q dependence of the transparency when using two 
different spectral functions and two different Fermi dis- 
tributions in the simulation, and the Q 2 dependent un- 



certainty from reactions mechanisms not included in the 
simulation. The uncertainty from the reaction mecha- 
nisms not included in the simulation were estimated by 
quantifying the difference in shape of the missing-mass 
spectra from data and simulation. 

Table ITVl lists the systematic uncertainties in extract- 
ing the unseparated cross sections and the uncertainty 
in extracting the ratio of longitudinal to transverse cross 
section. These were obtained using the quadrature sum 
of the correlated, uncorrelated and A dependent un- 
certainties along with two additional sources of uncer- 
tainty labeled as "Kinematics" and "momentum cover- 
age." The "Kinematics" part represents the uncertain- 
ties in the knowledge of spectrometer angles and momen- 
tum setting. It will not contribute to the nuclear trans- 
parency, since the spectrometer settings are exactly the 
same for different targets. The "momentum coverage" 
part represents the effect of the wider momentum cut 
(in order to obtain enough statistics in the low e region) 
used in the Rosenbluth separation than the extraction of 
nuclear transparency. 



IX. RESULTS AND DISCUSSION 

The differential cross sections for hydrogen are shown 
in Fig. [13] and those for all four nuclear targets (deu- 
terium, carbon, copper and gold) are shown in Fig. 1141 
All numerical values are tabulated in Appendix A and B. 
In the following sections, the global dependencies of the 
hydrogen and nuclear cross sections will be reviewed and 
the data compared to recent model calculations. The re- 
sults of the quasi-free reaction mechanism test and the 
extracted nuclear transparencies are also shown. 



A. Global dependencies and model comparison for 
hydrogen 

For the different values of Q 2 , the differential cross 
section shows the characteristic fall-off with —t, which 
may be due to the pion pole in <jl . The magnitude of the 
cross section decreases at constant W and with increasing 
Q 2 , mostly because the value of —t m i„ increases with Q 2 . 
In Fig. [T3] (bottom right panel) we show the differential 
cross section vs pion lab momentum (P w ). In this plot 
each point at a particular P^, represents the differential 
cross section averaged over the \t\ range shown in the one 
of the other panels of Fig. [12] that corresponds to the Q 2 
values at that P v (see Tabic [J). 

The cross sections are compared to predictions of 
two different models of pion electroproduction, the 
VGL- Regg e model [45| . and the more recent "KGM" 
model [46|. The VGL-Regge calculations are in a gauge 
invariant model incorporating 7r and p Rcggc trajectory 
exchanges. They significantly underestimate the mea- 
sured differential cross sections. Most of the discrep- 
ancy can likely be attributed to the model underesti- 
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FIG. 13: Differential cross sections % forpion electro-production from hydrogen versus \t\. The last panel (bottom right) shows 
the differential cross section versus the pion lab momentum (only the seven W =2.2 GeV points are shown here). For each of 
the points shown in the cross section versus pion momentum plot (last panel) the data were averaged over the respective t range 
shown in the previous panels. For the panels showing differential cross section versus \t\ the center of mass energy is W = 2.2 
GeV for all except one kinematic setting where W = 1.8 GeV (bottom middle). The inner error bars represent the statistical 
uncertainties, while the outer error bars are the sum in guadrature of the statistical and systematic uncertainties. The data are 
compared with both the VGL-Regge fldl l and the KGM JTdll calculations where available (see text). 



mating the transverse part of the cross section as shown 
in Ref. [47j], while they agree well with ol- The recent 
KGM model [46|, which includes a deep-inelastic scat- 
tering (DIS) Ansatz for the transverse part of the cross 
section, with the longitudinal cross section dominated by 
hadronic degrees of freedom and the pion electromag- 
netic form factor, agrees much better with the measured 
differential cross sections. 



B. Global dependencies for nuclear targets 

We have extracted the differential cross sections for 
all four nuclear targets (deuterium, carbon, copper and 
gold) at the eight different kinematics settings given in 
Table I. Here, the additional complication introduced due 
to the added degree of freedom induced by the Fermi- 
motion (or more general, the nuclear binding) of the 
struck proton is taken into account by extracting the 
doubly-differential cross sections dt &puu > where P^ M is 



the pion momentum in the center-of-mass frame of the 
virtual photon and the nucleus. 

The local variations in the nuclear cross sections as il- 
lustrated in Fig. [TJ] indicate effects due to Fermi motion 
beyond the central kinematics. These local variations 
are more pronounced for the deuterium target, because 
of its narrower Fermi cone. Although the general trend of 
the nuclear cross section is similar to that of the hydro- 
gen cross section, the fall-off of the nuclear cross sections 
with P v is steeper than that found for the hydrogen cross 
sections. 



C. Verification of the quasi- free mechanism 

A prerequisite for an interpretation of nuclear trans- 
parency as a function of Q 2 is that the reaction mech- 
anism remains identical over the Q 2 range. This trans- 
lates into an important condition in searching for CT us- 
ing pion electroproduction: the reaction should proceed 
through a quasi-free mechanism. 
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FIG. 14: T/ie extracted doubly- differential cross sections dtd d pCki versus the pion momentum P n for the deuterium, carbon, 

copper, and gold targets, respectively. The cross sections are normalized by Z, since 7r + can only be generated from a proton in 
this channel. The inner error bars represent the statistical uncertainties, while the outer error bars are the sum in quadrature of 
the statistical and systematic uncertainties. For each target at Q 2 =2.15 GeV 2 and Q 2 = 3.91 GeV 2 the solid symbol represents 
the high e kinematics while the open symbols represent the low e kinematics. 



There are many mechanisms that could break down 
the quasi-free assumption, such as: 

• Nucleon-nucleon (NN) final-state interac- 
tions The amplitude of NN final-state interactions 
can interfere with the elementary electro-pion pro- 
duction amplitude, which can change the ratio of 
the longitudinal and the transverse cross section. 
For example, the disagreement between the miss- 
ing mass spectra from data and simulation at Q 2 



= 1.1 (GeV) 2 (Figs Hand U is expected to be due 
to such effects, and other potential reaction mecha- 
nism effects not included in the Monte Carlo model. 
The agreement between data and Monte Carlo at 
higher values of Q 2 , however, suggests only a small 
contribution from such NN final-state interaction 
complications. 

• Rescattering Rcscattering involves the electro- 
production of a meson followed by a second inter- 
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action that produces the detected ir + particle. For 
example, such rescattering contributions have been 
shown to dominate the cross section in p° photo- 
production for t > 0.5 (GeV) 2 [isf . In principle, 
if the rescattering effect dominates, one would ex- 
pect a modification of both the longitudinal to the 
transverse cross section ratio, and the cross section 
dependence on W. 

• Pion excess Excess pions may be present in a 
nuclear system due to the long range of meson- 
exchange currents [50j . If such pion excess effects 
are significant, one would anticipate a change in the 
ratio ctl/ctt of the measured nuclear cross sections 
as compered to the hydrogen cross sections. An 
earlier experiment on light nuclei at low Q 2 did not 
find any pion excess [49| . 

• Medium modification of nucleons The 

European-Muon Collaboration (EMC) discovered 
[5~fl ] that the structure functions for deep inelastic 
inclusive lepton scattering off nuclear targets 
differed from deutcron targets. Although the effect 
remains poorly understood, it is generally accepted 
that nuclcon structure will be modified within a 
nuclear medium, both due to nuclear binding and 
to non-nuclconic QCD effects. Such a medium 
modification of nucleons could also impact the 
nuclear pion electro-production cross sections, 
but is generally expected to also lead to a change 
in the ratio of longitudinal and transverse pion 
electro-production cross sections, (Jl/ctt- 

• Two-nucleon correlations A series of A(e, e'p) 
measurements revealed that spectroscopic factors 
for proton valence shells were quenched by approx- 
imately 30-35% as compared to mean-field expec- 
tations [13, H3|- A possible explanation for this 
discrepancy is that correlations move some of the 
single-particle strength to orbitals above the Fermi 
energy. This kind of correlation will change the 
nucleus spectral function, and thus break down the 
quasi-free assumption. 

The most straightforward verification of the quasi- 
free mechanism is the equivalence of the longitudinal- 
transverse character for pion electroproduction from nu- 
clear and hydrogen targets. The ratio of longitudinal to 
transverse cross sections should, for instance, be inde- 
pendent of the nuclear atomic number A. The ratio of 
longitudinal to transverse cross sections at fixed Q 2 = 
2.15 and 3.91 GeV 2 for the various targets in this exper- 
iment are shown in Fig. 1121 

The A-dependence of the ctl/&t ratio agrees with a 
flat line fit within the experimental uncertainties, and are 
thus consistent with the quasi-free assumption. However, 
they cannot rule out non quasi-free reaction mechanism 
effects that affect the longitudinal and transverse char- 
acter of pion electroproduction in a similar fashion, for 
instance, NN final-state interaction, rescattering, pion 



excess, or medium modification effects. Note that we 
have intentionally kept the value of —t of the measure- 
ments low (< 0.5 GeV) to minimize complications due 
to rescattering or two-nucleon effects. Together with the 
overall good agreement of data and Monte Carlo simula- 
tions (beyond Q 2 = 1.1 GeV 2 ), we have gained confidence 
in the validity of the quasi-free reaction mechanism. 

D. Nuclear Transparencies 

As mentioned earlier, nuclear transparency is defined 
as the ratio of the cross section per nuclcon for a process 
on a bound nucleon inside a nucleus to that from a free 
nucleon. The P w dependence of the nuclear transparen- 
cies is shown in Fig. 1151 Results are shown for deuterium 
and carbon (top left), copper (bottom left) and gold (bot- 
tom right) and aluminum (top right), which was also used 
for determining the cell wall background contribution for 
cryogenic targets. 

Before we discuss the various model calculations, we 
first redefine the nuclear transparency (To) as the cross 
sections of heavy nuclear targets as compared to those 
of a deuterium target. This reduces the uncertainty due 
to both the unknown elementary pionclcctroproduction 
off a neutron and effects due to the procedure to take 
into account Fermi motion. The results are shown in 
Fig. [lHl where we present the data versus Q 2 rather than 
pion momentum P„.. Of course, the results are not too 
different, as could be anticipated from the nuclear trans- 
parency results of deuterium to be close to unity, as can 
be seen in top left panel of Fig. Qj)] We note in addi- 
tion that the deuteron nuclear transparency results are 
found to be consistent with a constant-line fit with 81% 
probability. The nuclear transparency results are listed 
in Appendix C. 

The nuclear transparencies are expected to be near- 
constants over the pion momentum range of the experi- 
ment from traditional nuclear physics point of view [54| , 
because the hadron-nucleon cross sections are nearly in- 
dependent of momentum over the range of momenta in 
Figs.fT51and[TBl Instead, the observed pion nuclear trans- 
parency results (as compared both to hydrogen and deu- 
terium cross sections) shows a slow but steady rise versus 
pion momentum for the nuclear (A > 2) targets, causing 
a deviation from calculations without CT. 



E. Comparison with model calculations 

We compare our results with the calculations of Lar- 
son, et al. [f| (solid and dashed curves) , Cosyn, et al. (55| 
(dot-short dash and dot-long dash curves) and Kaskulov, 
et al. 56] (dotted and dot-dot-dash). 

Larson, et al. compute the nuclear transparency at 
the exact kinematics of the experiment in terms of semi- 
classical formula based on the eikonal approximation and 
a parametrization of the effects of final state interactions 
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FIG. 15: (color online) Nuclear transparency, T, vs P w for 
2 H and 12 C, 27 Al, 63 Cu, and 197 Au. The inner error bars are 
the statistical uncertainties and the outer error bars are the 
statistical and point-to-point systematic uncertainties added 
in quadrature. The solid circles (blue) are the high e points 
while the solid squares (red) are low e points. The dark (blue) 
bands are the model uncertainties (for details see Sec. I VW\) . 
The dashed and solid lines (red) are Glauber calculations from 
Larson, et al. Jj^j, with and without CT respectively. Similarly, 
the dot-short dash and dot-long dash lines (blue) are Glauber- 
calculations with and without CT from Cosyn, et al. JEh l. The 
effects of short-range correlations are included in these latter 
calculations. The dotted and dot- dot- dashed lines (green) are 
microscopic-h BUU transport calculations from Kaskulov et 
al. fEdl l. with and without CT respectively. 



(FSI) in terms of an effective interaction. This semi- 
classical formula involves a single integral over the path of 
the outgoing pion which is suited for situations in which 
the kinematics of the final pion are known. The nuclear 
density is taken as a Woods-Saxon form with radius pa- 
rameter R = 1.1 A 1 / 3, fm and diffuscncss a = 0.54 fm. 
The effective interaction is based on the quantum diffu- 
sion model of Ref. 0], which predicts the interaction of 
the PLC to be approximately proportional to the prop- 
agation distance z for z < l c The coherence length (or 
formation length) is parametrized as described earlier in 
section [TTJ Larson, Miller and Strikman use the following 
parameters; Tiif etime = 1 fm/c and M\ = 0.7 (GeV/c 2 ) 2 . 
In the limit of l c = a PLC is not created and the ef- 
fective interaction reduces to a Glauber-type calculation 
with (J e ff « ovaK-FV), the 7r — N cross section for pion 
momentum P„ obtained from a parametrization by the 
Particle Data Group (PDG) [57|. 

On the other hand, Cosyn et al. calculate the nu- 



clear transparency as a ratio of the eight fold differen- 
tial cross section for pion electroproduction in a relativis- 
tic multiple-scattering Glauber approximation (RMSGA) 
integrated over the kinematic range of the experiment 
to that in a relativistic plane wave impulse approxima- 
tion (RPWIA). In the RPWIA all particles are taken to 
be relativistic plane waves, while in the RMSGA, the 
wave function of the spectator nucleon and the outgo- 
ing pion is taken to be a convolution of a relativistic 
plane wave and a Glauber-type eikonal phase operator 
that parametrizes the effects of FSI. The parametriza- 
tion chosen by these authors reflect the diffractive nature 
of nucleon-nucleon (N'N) and pion-nucleon (wN) colli- 
sions at intermediate energies. The parameters er*^ (to- 
tal cross section, with i as the outgoing pion or nucleon), 
PiN (slope parameter), and e^jv (ratio of real to imagi- 
nary part of the scattering amplitude) were obtaine d by 
fits to the N'N -> N'N databases from the PDG [57[, 
for i = N'. For i = n, the parameters were obtained 
from fits to PDG [57| databases, SAID [HI and Ref. [59[. 
For outgoing nuclcons with kinetic energy lower than ~ 
300 MeV, the Glauber formalism is no longer applicable 
and the FSI were parametrized in a relativistic optical 
model eikonal approximation (ROMEA) [6Cj], with the 
global (S-V) optical model parametrization of Cooper et 
al. [6l| . CT was incorporated by replacing the total cross 
section parameterc*^ with an effective one based on the 
quantum diffusion model this mirrors the effective 
interaction parameter of Larson et al., described earlier. 
The parameters used for l c were exactly the same for 
both set of authors. Cosyn et al. also include the effects 
of short range correlations (SRC) in their calculations. 
The Glauber phase factor described above is corrected 
for SRC by replacing the single nucleon density typically 
used in Glauber-type calculations with an effective den- 
sity. The effective density modifies the single nucleon 
density with a Jastrow correlation function and normal- 
ization functions which ensure the integral of the effective 
density is equal to the total number of nucleons. 

We have also compared our results to the recent cal- 
culations of Kaskulov et al., who calculate nuclear trans- 
parency as the ratio of their model differential cross sec- 
tion calculated in the lab frame, with and without FSI, 
where both types of model cross sections are integrated 
over the kinematic range of the experiment. Their model 
is built around a microscopic description [46] of the ele- 
mentary process p(e, e'ir + )n which is divided into a soft 
hadronic and a hard partonic or deep inelastic scatter- 
ing production. For the reaction on nuclei, the ele- 
mentary interaction is kept the same, and standard ef- 
fects such as Fermi motion, Pauli blocking and nuclear 
shadowing are accounted for. Finally all produced pre- 
hadrons and hadrons are propagated through the nuclear 
medium according to the Boltzmann-Uchling-Uhlcnbcck 
(BUU) transport equation. The DIS contribution to the 
cross section is determined by the Lund fragmentation 
model [62j and the time development of the interactions 
of the pre-hadron is determined by the quantum diffusion 
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FIG. 16: (color online) The redefined Nuclear transparency 
T D versus Q 2 for 12 C, 27 Al, 63 Cu, and 197 Au. Curves rep- 
resent calculations identical to those in Fig. \15\ Model un- 
certainties are identical to Fig. \15\ are hence are not repeated 
here. 



model [4|. The production time and the formation time 
are from a Monte Carlo calculation based on the Lund 
fragmentation model [62| described in Ref. (63[. Only 
the DIS part of the cross section is effected by this pre- 
hadronic interaction and thus in this model only the DIS 
events are responsible for the CT effect. 

Our results are in good agreement with the CT cal- 
culations of Larson et al., while the calculations of both 
Cosyn et al. and Kaskulov et al. seem to somewhat 
overestimate the data. However, whereas one can argue 
about details of the calculation, it is more important to 
note that the trend of all calculations including CT in 
Figs. IT51 and [TB]are consistent with the trend of the nu- 
clear transparency data versus Q 2 . 

The underlying cause for the rise in nuclear trans- 
parency is different for the different model calculations, 
however. Whereas the dominance of the longitudinal- 
photon production mechanism in exclusive (e,e'7r + ) reac- 
tion is thought to proceed only at asymptotic Q 2 , and one 
would thus anticipate to see the CT effects enter this lon- 
gitudinal channel, the Kaskulov et al. calculations find 
the CT effects coming from the transverse-photon pro- 
duction mechanism. Thus, we separately show in Fig. 1151 
the nuclear transparency results for our low e data (but 
note that the calculations shown in Fig. [TB] are for the 
high e kinematics only). Clearly, within the experimen- 
tal uncertainties, we do not see any obvious difference 
between the nuclear transparency results measured at 
high and at low e values. Given the uncertainties of the 



present experiment, we can not distinguish between the 
suggested mechanisms. 



F. Further studies of the CT mechanism 

Figs. [15] and [16] show a rise of nuclear transparency 
with outgoing pion momentum, or alternatively Q 2 , 
which deviates from the traditional nuclear physics ex- 
pectation. However, as can be readily witnessed from 
Table [I] there exists a strong correlation between the 
outgoing pion momentum, P n , and the magnitude of the 
virtual-pion (three-)momentum, k^. This poses a poten- 
tial pitfall in that the observed CT-like behavior could 
be an artifact due to increased values of k n (and thus in- 
creased probability for reaction mechanisms beyond the 
quasi- free picture), rather than a dependence on P n . To 
investigate this further, we performed measurements at 
different values of k n (and thus different P v ), but at iden- 
tical Q 2 (2.15 GeV 2 ). 

The results are shown in Fig. [T7l and indicate that the 
nuclear transparency does not show any obvious depen- 
dence on k^. This result rules out the possibility that 
nuclear transparency only depends on k v (but does not 
yet rule out the possibility that nuclear transparency de- 
pends on k n as well as other variables). 

The values of Q 2 and P n for this exclusive (e,e'7r + ) 
experiment are also strongly correlated. Since Q 2 is in 
general terms related to the size of the PLC, and P w 
to the formation length of the PLC, the rise in nuclear 
transparency results can also be a mixed PLC size and 
formation-length effect. To further study which of these 
two effects dominates, we formed nuclear transparency 
ratios of the heavy target nuclei ( 27 A1, 63 Cu, and 197 Au) 
with respect to 12 C termed as Tc- This ratios should 
be less sensitive to formation length effects. Here we use 
the nuclear size as a yardstick to gauge against formation 
length effects. 

The "super ratios" T c for 27 Al (top right panel), 63 Cu 
(bottom left panel), and 197 Au (bottom right panel) are 
shown in Fig. [T8l The results are consistent with a flat 
line, within the (large) experimental uncertainties, with 
probabilities of 0.32, 0.40, and 0.64, respectively. A plau- 
sible explanation is that the pion formation length in our 
kinematics is already much larger than the nucleus ra- 
dius (a simple estimate would give a formation length of 
the order of 10 fm). We also note that a reasonable ap- 
proximation for the coherence length of the virtual pho- 
ton (l c = l/(2Mx)) renders a value of about 0.2-0.5 fm, 
which is much smaller than the size of the nucleus. Thus, 
we should not be sensitive to any coherence length effects. 
We conclude that these results seem to favor the expla- 
nation that the rise in nuclear transparency measured is 
due to a small PLC size, rather than a PLC formation 
length effect. 
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FIG. 17: (color online) Nuclear transparency, T vs k-n for 2 H, 
12 C and 27 Al. 
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G. The A dependence 

The dependence of the nuclear transparency data on 
the atomic number A gives further insight in the proper 
interpretation of the data in terms of an onset of CT. This 
goes beyond the Q 2 or P x dependence of nuclear trans- 
parencies for one single nuclear target described above. 
Here, the entire nuclear transparency data set was exam- 
ined using a simple description with an effective param- 
eter a (for each value of Q 2 ). Using only one parameter 
obviously neglects specific surface effects of the various 
nuclei, but has proven to be an effective way to describe 
bulk properties of the nuclear medium. For example, 
pion-nucleus scattering total cross section data are well 
described using such a single parameter, as a A = A a a N , 
with a the nuclear cross section, a N the nucleon cross 



Q (GeV ) 



section, and a = 0.76 [64[. No noticeable dependence on 
the incident pion energy was measured. 

The full nuclear transparency data set from pionCT 
was fitted as function of Q 2 (with T = A a_1 which 
equates to the same fit form as above) are shown in 
Fig. 1191 The uncertainties are dominated by systematics, 
and include contributions from both fitting and model 
uncertainties. The results indicate a value for the pa- 
rameter a clearly deviating from the total pion-nucleus 
scattering cross section results, with values a > 0.76. 
Furthermore, a noticeable dependence of a on Q 2 , or 
equivalently the pion kinetic energy, was measured. 

We compare the extracted values of a, as function of 
Q 2 , with the calculations including CT effects of Larson 
et al. @ and Cosyn et al. (H- The agreement with the 
calculations from Ref. || is excellent, but the data are 
systematically below the calculations (including both CT 
and short-range correlation effects) of Ref. [55[. As men- 
tioned above, Kaskulov et al. [46| recently suggested that 
the CT effect should only exist in the transverse cross sec- 
tion. For this reason, we have again separately indicated 
the a values for the low and high e values of this exper- 
iment, where applicable. Within our uncertainties, we 
see no indication of this prediction. However, we have 
to warn that the difference between the low- to high-e 
kinematics in terms of contributions from the longitudi- 
nal cross sections to the total measured cross sections 
changes by less than 30% in our kinematics (see Fig. [T^] 
and Table It is thus very well possible that the mea- 
sured effect is solely due to the transverse contributions. 



FIG. 18: (color online) The Tc vs Q 2 , where Tc is defined 
as the super ratio of the heavy nuclear targets with respect 
to 12 C (T(A > 12) to T(A = 12)). The inner error bars 
are the statistical uncertainties, the outer error bars are the 
statistical and point-to-point systematic uncertainties added in 
quadrature. The probabilities of the fitting to a straight line 
are 32%, 40% and 64% for carbon, copper and gold targets, 
respectively. The calculations are same as in Fig. \15\ but 
divided by the results for 12 C. 



X. CONCLUSIONS 

The A(e, e'7r + ) reaction was studied on a variety of nu- 
clear targets: 1 H, 2 H, 12 C, 27 Al, Cu, and Au. Data were 
taken up to a four-momentum transfer squared of Q = 
4.8 GeV 2 , and analyzed in terms of nuclear transparen- 
cies, the escape probability of the positively-charged pion 
from the nuclear medium. A rise of the nuclear trans- 
parency with Q 2 or pion momentum could signal an on- 
set of Color Transparency, expected to occur at large val- 
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FIG. 19: (color online) The parameter a, as extracted from 
the global nuclear data set of this experiment (from T = 
versus Q (solid black circles). The inner error bars 
indicate the statistical uncertainty, and the outer error bars 
are the quadrature sum of statistical, systematic and model- 
ing uncertainties. The hatched band is th e value of a extracted 
from pion-nucleus scattering data '6Aj. The solid, dashed 
and dotted lines are a obtained from fitting the A-dependence 
of the theoretical calculations: the Glauber and Glauber+CT 
calculations of Ref. f^j, and the Glauber + CT (including 
short-range correlation effects) calculations of Ref. HI/, re- 
spectively. The red solid squares in addition show the a value 
extracted at low e value at two values of Q 2 . 



ucs of Q 2 from both pcrturbativc and non-perturbativc 
QCD. 

The existence of Color Transparency is well-known. 
The most convincing evidence is the analysis of Fermi- 
lab E791 data on the diffractive dissociation in two jets. 
However, conclusive proof of the onset of CT is more elu- 
sive, although hints of it were recently seen in the analysis 
of p vector meson data at fixed coherence length [2l[ and 
7r~ photo-production [23| . Our results show a far more 
conclusive onset of CT over the relatively large range in 
pion momentum between 2.5 and 4.5 GeV/c, and the 
Q 2 range between 1.1 and 4.8 GeV 2 , and are in good 
agreement with model calculations including the CT ef- 
fect of Larson et al., Cosyn et al., and Kaskulov et al. 
0, HH, [56| . The latter calculation also agrees excellently 
with the measured —t dependence of the differential cross 
sections, giving further credence to the onset of CT. 

Specialized data sets were added to ensure that the 
noted rise of nuclear transparencies is indeed due to CT. 
In particular, the cross checks performed are: 

1. L/T character of the cross section 

Within uncertainties, we find the longitudinal- 
transverse character of the pion electroproduction 
cross sections to be similar off protons and heavier 
nuclei. This supports a quasifree reaction mecha- 
nism. 



2. Virtual-pion momentum k„ 

The measured nuclear transparencies do not appear 
to depend on the virtual-pion (three-)momentum, 
which could be a signal for increased reaction mech- 
anism effects beyond the quasifree picture. This 
suggests that reaction mechanism effects beyond 
the quasifree are suppressed. 

3. Cross section ratios of medium- to heavy- nuclei 
"Super ratios" of the heavier target nuclei ( 27 A1, 
Cu, Au) with respect to 12 C show a similar rise 
in nuclear transparencies, indicating that the pion 
formation length in our kinematics is already far 
longer than the nuclear radius. 

4. A dependence 

The A dependence of the data can be described by a 
single parameter, a. It was found to be consistently 
larger than the 0.76 found from pion-nucleus total 
cross section data, and rising with Q 2 , consistent 
with a CT Ansatz. 

Furthermore, our results are at nearly-constant and small 
coherence length (« 0.4 fm), such that possible compli- 
cations due to irp exchange terms are minimized. 

Our results are consistent with the predicted early on- 
set of CT in mesons compared to baryons, and together 
with previous meson transparency measurements (2ll . [23[ 
suggest a gradual transition to meson production with 
small inter-quark separation. These results put severe 
constraints on early models of CT which predict a dra- 
matic transition with a threshold-like behavior. The un- 
ambiguous observation of the onset of CT uniquely points 
to the role of color in exclusive high-Q 2 processes. Fur- 
thermore, it is an effective signature of the approach to 
the factorization regime in meson electroproduction ex- 
periments, necessary for the access to Generalized Parton 
Distributions. These results will be further extended by 
data to be taken after the Jefferson Lab upgrade to 12 
GeV, with planned exclusive A(e, e'ir + ) measurements up 
to Q 2 « 10 GeV 2 [H. 
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Appendix A 



Target 


Q 2 


W 


-t 


Pit 


e 


pUM 


da 


stat. err. 


sys. err. 




GeV 2 


GeV 


GeV 2 


GeV/c 




GeV/c 


fib/GeV* 


fib/GeV' s 


Hb/GeV' 6 


LH 2 


1.71 


2.36 


0.100 


3.35 


0.532 


0.99 


3.72 


0.11 


0.24 


LH 2 


1.80 


2.33 


0.117 


3.31 


0.538 


0.97 


3.46 


0.11 


0.23 


LH 2 


1.84 


2.32 


0.128 


3.30 


0.539 


0.96 


3.47 


0.11 


0.20 


LH 2 


1.87 


2.31 


0.136 


3.30 


0.538 


0.96 


3.04 


0.10 


0.19 


LH 2 


1.90 


2.30 


0.144 


3.26 


0.545 


0.95 


3.07 


0.10 


0.19 


LH 2 


1.93 


2.28 


0.152 


3.24 


0.548 


0.94 


2.92 


0.09 


0.17 


LH 2 


1.94 


2.27 


0.160 


3.21 


0.554 


0.93 


2.96 


0.09 


0.18 


LH 2 


1.95 


2.27 


0.167 


3.21 


0.555 


0.93 


2.78 


0.09 


0.17 


LH 2 


2.02 


2.25 


0.174 


3.20 


0.554 


0.92 


2.75 


0.08 


0.16 


LH 2 


2.06 


2.24 


0.181 


3.19 


0.553 


0.91 


2.54 


0.08 


0.15 


LH 2 


2.05 


2.24 


0.187 


3.19 


0.553 


0.91 


2.71 


0.09 


0.15 


LH 2 


2.09 


2.23 


0.193 


3.17 


0.555 


0.90 


2.61 


0.09 


0.14 


LH 2 


2.08 


2.24 


0.200 


3.17 


0.555 


0.91 


2.48 


0.07 


0.13 


LH 2 


2.09 


2.23 


0.208 


3.17 


0.554 


0.90 


2.38 


0.07 


0.12 


LH 2 


2.13 


2.23 


0.216 


3.16 


0.552 


0.90 


2.55 


0.07 


0.13 


LH 2 


2.13 


2.23 


0.224 


3.16 


0.552 


0.90 


2.27 


0.07 


0.11 


LH 2 


2.13 


2.23 


0.232 


3.16 


0.552 


0.90 


2.29 


0.07 


0.11 


LH 2 


2.20 


2.19 


0.240 


3.10 


0.563 


0.88 


2.21 


0.07 


0.11 


LH 2 


2.22 


2.19 


0.249 


3.10 


0.560 


0.87 


2.18 


0.06 


0.11 


LH 2 


2.19 


2.19 


0.259 


3.10 


0.563 


0.88 


2.15 


0.07 


0.11 


LH 2 


2.26 


2.18 


0.270 


3.08 


0.558 


0.86 


2.11 


0.06 


0.11 


LH 2 


2.29 


2.17 


0.283 


3.08 


0.557 


0.86 


2.05 


0.06 


0.12 


LH 2 


2.34 


2.16 


0.299 


3.07 


0.557 


0.86 


1.95 


0.06 


0.12 


LH 2 


2.44 


2.13 


0.317 


3.06 


0.556 


0.84 


2.11 


0.07 


0.15 


LH 2 


2.53 


2.12 


0.341 


3.04 


0.547 


0.82 


1.78 


0.06 


0.04 


LH 2 


2.63 


2.08 


0.385 


3.03 


0.561 


0.82 


1.83 


0.07 


0.19 


LH 2 


3.75 


2.21 


0.425 


3.89 


0.251 


0.89 


0.64 


0.02 


0.06 


LH 2 


4.00 


2.14 


0.548 


3.80 


0.255 


0.85 


0.58 


0.02 


0.06 


LH 2 


4.00 


2.14 


0.646 


3.73 


0.255 


0.84 


0.64 


0.02 


0.07 


LH 2 


4.21 


2.07 


0.758 


3.63 


0.261 


0.80 


0.70 


0.02 


0.07 


LH 2 


2.55 


2.28 


0.199 


3.55 


0.435 


0.94 


1.82 


0.06 


0.09 


LH 2 


2.64 


2.26 


0.236 


3.52 


0.433 


0.93 


1.70 


0.05 


0.09 


LH 2 


2.68 


2.24 


0.260 


3.49 


0.438 


0.91 


1.50 


0.05 


0.08 


LH 2 


2.72 


2.23 


0.280 


3.46 


0.441 


0.91 


1.41 


0.05 


0.07 


LH 2 


2.78 


2.21 


0.300 


3.45 


0.440 


0.90 


1.43 


0.05 


0.08 


LH 2 


2.85 


2.19 


0.320 


3.42 


0.443 


0.88 


1.42 


0.04 


0.08 


LH 2 


2.90 


2.18 


0.340 


3.40 


0.445 


0.87 


1.26 


0.04 


0.07 


LH 2 


2.94 


2.16 


0.360 


3.38 


0.445 


0.86 


1.36 


0.04 


0.07 


LH 2 


2.96 


2.16 


0.380 


3.37 


0.446 


0.86 


1.28 


0.04 


0.07 


LH 2 


3.06 


2.13 


0.400 


3.33 


0.450 


0.84 


1.26 


0.04 


0.07 


LH 2 


3.08 


2.12 


0.420 


3.30 


0.454 


0.83 


1.22 


0.04 


0.06 


LH 2 


3.11 


2.11 


0.445 


3.29 


0.454 


0.82 


1.16 


0.03 


0.06 


LH 2 


3.27 


2.07 


0.474 


3.28 


0.450 


0.80 


1.08 


0.04 


0.06 


LH 2 


3.25 


2.08 


0.508 


3.25 


0.447 


0.80 


1.09 


0.04 


0.06 


LH 2 


3.50 


1.99 


0.565 


3.21 


0.455 


0.77 


1.22 


0.04 


0.07 


LH 2 


0.92 


2.31 


0.038 


2.84 


0.490 


0.96 


10.74 


0.32 


0.53 


LH 2 


0.98 


2.29 


0.047 


2.81 


0.497 


0.94 


9.56 


0.26 


0.49 


LH 2 


1.00 


2.28 


0.053 


2.80 


0.499 


0.94 


8.99 


0.26 


0.46 



TABLE V: Extracted cross sections and their uncertainties for hydrogen target at eight kinematics settings. The Q 2 , W, —t, 
P-n and e values for each kinematics are shown at the same time. 



Target 


Q 


W 


-t 




e 


pCM 


da 


stat. err. 


sys. err. 




GeV^ 


GeV 


GeV^ 


GeV/c 




GeV/c 


fib/GeV 


fj,b /GeV' s 


[ib 1 GeV 6 


LH 2 


1.04 


2.27 


0.058 


2.79 


0.501 


0.93 


8.15 


0.22 


0.42 


LH 2 


1.03 


2.27 


0.062 


2.78 


0.501 


0.93 


7.87 


0.23 


0.42 


LH 2 


1.10 


2.26 


0.066 


2.79 


0.497 


0.92 


7.90 


0.22 


0.42 


LH 2 


1.07 


2.26 


0.070 


2.76 


0.504 


0.92 


7.35 


0.21 


0.39 


LH 2 


1.10 


2.25 


0.074 


2.77 


0.502 


0.92 


6.73 


0.19 


0.35 


LH 2 


1.10 


2.25 


0.078 


2.76 


0.501 


0.92 


7.06 


0.19 


0.37 


LH 2 


1.12 


2.24 


0.082 


2.75 


0.503 


0.91 


6.71 


0.19 


0.34 


LH 2 


1.13 


2.24 


0.086 


2.74 


0.506 


0.91 


6.52 


0.19 


0.33 


LH 2 


1.16 


2.23 


0.090 


2.74 


0.504 


0.91 


6.34 


0.17 


0.32 


LH 2 


1.16 


2.24 


0.095 


2.73 


0.503 


0.91 


6.00 


0.17 


0.30 


LH 2 


1.19 


2.22 


0.100 


2.71 


0.510 


0.89 


6.58 


0.19 


0.33 


LH 2 


1.20 


2.22 


0.105 


2.72 


0.505 


0.90 


6.01 


0.19 


0.31 


LH 2 


1.21 


2.22 


0.110 


2.71 


0.505 


0.89 


6.50 


0.20 


0.33 


LH 2 


1.18 


2.23 


0.116 


2.70 


0.506 


0.90 


6.18 


0.19 


0.32 


LH 2 


1.16 


2.24 


0.123 


2.71 


0.501 


0.90 


5.51 


0.18 


0.30 


LH 2 


1.17 


2.23 


0.133 


2.69 


0.505 


0.89 


6.15 


0.18 


0.37 


LH 2 


1.21 


2.22 


0.154 


2.67 


0.502 


0.88 


7.21 


0.20 


0.59 


LH 2 


2.06 


2.24 


0.165 


3.21 


0.268 


0.91 


2.52 


0.08 


0.12 


LH 2 


2.16 


2.21 


0.209 


3.17 


0.269 


0.90 


2.01 


0.07 


0.10 


LH 2 


2.26 


2.18 


0.248 


3.12 


0.274 


0.88 


1.95 


0.06 


0.10 


LH 2 


2.27 


2.18 


0.290 


3.07 


0.276 


0.86 


1.97 


0.06 


0.10 


LH 2 


2.28 


2.18 


0.326 


3.02 


0.274 


0.85 


2.04 


0.07 


0.13 


LH 2 


1.81 


1.86 


0.241 


2.21 


0.634 


0.68 


5.74 


0.19 


0.43 


LH 2 


1.86 


1.84 


0.276 


2.17 


0.634 


0.67 


5.01 


0.16 


0.33 


LH 2 


1.93 


1.82 


0.305 


2.15 


0.633 


0.65 


5.10 


0.15 


0.31 


LH 2 


2.00 


1.79 


0.330 


2.12 


0.636 


0.63 


5.20 


0.18 


0.31 


LH 2 


2.08 


1.76 


0.355 


2.09 


0.638 


0.62 


5.59 


0.15 


0.32 


LH 2 


2.14 


1.74 


0.385 


2.07 


0.637 


0.60 


5.61 


0.15 


0.32 


LH 2 


2.19 


1.74 


0.414 


2.07 


0.629 


0.60 


5.73 


0.16 


0.33 


LH 2 


2.25 


1.71 


0.444 


2.04 


0.630 


0.58 


5.58 


0.18 


0.30 


LH 2 


2.30 


1.70 


0.479 


2.02 


0.628 


0.57 


5.61 


0.19 


0.30 


LH 2 


2.40 


1.67 


0.548 


1.98 


0.622 


0.55 


5.58 


0.17 


0.29 


LH 2 


4.39 


2.32 


0.469 


4.49 


0.259 


0.96 


0.45 


0.02 


0.02 


LH 2 


4.71 


2.24 


0.621 


4.37 


0.265 


0.91 


0.41 


0.01 


0.02 


LH 2 


4.82 


2.21 


0.765 


4.28 


0.269 


0.90 


0.41 


0.01 


0.02 


LH 2 


4.98 


2.17 


0.914 


4.17 


0.267 


0.86 


0.43 


0.01 


0.02 


LH 2 


3.45 


2.38 


0.280 


4.24 


0.381 


1.00 


0.75 


0.02 


0.04 


LH 2 


3.55 


2.36 


0.340 


4.20 


0.381 


0.99 


0.71 


0.02 


0.04 


LH 2 


3.69 


2.32 


0.380 


4.17 


0.383 


0.97 


0.67 


0.02 


0.04 


LH 2 


3.78 


2.30 


0.420 


4.12 


0.386 


0.95 


0.68 


0.02 


0.04 


LH 2 


3.84 


2.28 


0.460 


4.09 


0.386 


0.93 


0.62 


0.02 


0.03 


LH 2 


3.91 


2.26 


0.499 


4.05 


0.388 


0.92 


0.59 


0.02 


0.03 


LH 2 


3.92 


2.25 


0.540 


4.00 


0.393 


0.91 


0.62 


0.02 


0.03 


LH 2 


4.00 


2.23 


0.579 


3.97 


0.395 


0.90 


0.63 


0.02 


0.03 


LH 2 


4.10 


2.21 


0.619 


3.93 


0.395 


0.88 


0.61 


0.02 


0.03 


LH 2 


4.30 


2.15 


0.663 


3.91 


0.394 


0.85 


0.66 


0.02 


0.03 


LH 2 


4.40 


2.12 


0.741 


3.87 


0.398 


0.84 


0.63 


0.02 


0.03 



TABLE VI: Continue of Table. [V] 
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Appendix B 



Target 


Q 2 


W 


-t 


P-K 


pUM 


e 


da 
dtdP^ M 


stat. err. 


sys. err. 




GeV 2 


GeV 


GeV 2 


GeV/c 


GeV/c 




fib/GeV* 


fib/GeV' s 


tib/GeV' 6 


LD 2 


2.24 


3.68 


0.177 


3.20 


1.35 


0.559 


17.36 


0.63 


1.30 


LD 2 


2.24 


3.68 


0.188 


3.19 


1.34 


0.559 


31.72 


1.16 


1.82 


LD 2 


2.24 


3.68 


0.195 


3.18 


1.34 


0.559 


45.58 


1.73 


2.60 


LD 2 


2.24 


3.68 


0.203 


3.17 


1.34 


0.559 


55.31 


2.03 


3.10 


LD 2 


2.24 


3.68 


0.211 


3.16 


1.33 


0.559 


64.42 


2.48 


3.72 


LD 2 


2.24 


3.68 


0.220 


3.15 


1.33 


0.559 


66.01 


2.42 


3.98 


LD 2 


2.24 


3.68 


0.229 


3.14 


1.32 


0.559 


59.77 


2.20 


3.78 


LD 2 


2.24 


3.68 


0.240 


3.13 


1.32 


0.559 


53.28 


2.00 


3.51 


LD 2 


2.24 


3.68 


0.251 


3.12 


1.31 


0.559 


41.23 


1.55 


2.78 


LD 2 


2.24 


3.68 


0.263 


3.10 


1.31 


0.559 


27.73 


1.03 


1.86 


LD 2 


2.24 


3.68 


0.277 


3.08 


1.30 


0.559 


15.29 


0.61 


1.01 


LD 2 


2.24 


3.68 


0.299 


3.04 


1.28 


0.559 


7.86 


0.28 


0.46 


LD 2 


4.00 


3.86 


0.442 


3.90 


1.43 


0.254 


4.65 


0.21 


0.58 


LD 2 


4.00 


3.86 


0.512 


3.84 


1.41 


0.254 


10.20 


0.43 


1.25 


LD 2 


3.02 


3.73 


0.290 


3.46 


1.37 


0.449 


6.75 


0.25 


0.47 


LD 2 


3.02 


3.73 


0.316 


3.42 


1.36 


0.449 


17.55 


0.66 


1.13 


LD 2 


3.02 


3.73 


0.337 


3.40 


1.35 


0.449 


24.30 


0.92 


1.49 


LD 2 


3.02 


3.73 


0.359 


3.38 


1.34 


0.449 


25.32 


0.98 


1.50 


LD 2 


3.02 


3.73 


0.380 


3.36 


1.33 


0.449 


22.72 


1.00 


1.28 


LD 2 


3.02 


3.73 


0.402 


3.34 


1.32 


0.449 


16.21 


0.65 


0.91 


LD 2 


1.12 


3.60 


0.067 


2.78 


1.30 


0.504 


188.27 


6.96 


11.29 


LD 2 


1.12 


3.60 


0.071 


2.78 


1.30 


0.504 


262.48 


9.75 


14.97 


LD 2 


1.12 


3.60 


0.077 


2.77 


1.29 


0.504 


310.24 


11.76 


18.10 


LD 2 


1.12 


3.60 


0.082 


2.76 


1.29 


0.504 


316.58 


12.18 


19.00 


LD 2 


1.12 


3.60 


0.088 


2.75 


1.29 


0.504 


268.56 


10.35 


16.58 


LD 2 


1.12 


3.60 


0.094 


2.75 


1.29 


0.504 


211.69 


7.88 


13.73 


LD 2 


1.12 


3.60 


0.099 


2.74 


1.28 


0.504 


144.49 


5.42 


9.61 


LD 2 


1.12 


3.60 


0.105 


2.72 


1.27 


0.504 


75.80 


3.03 


4.91 


LD 2 


1.12 


3.60 


0.111 


2.70 


1.27 


0.504 


36.38 


1.60 


2.26 


LD 2 


1.12 


3.60 


0.120 


2.68 


1.25 


0.504 


16.95 


0.98 


0.95 


LD 2 


1.12 


3.60 


0.128 


2.66 


1.24 


0.504 


9.92 


0.69 


0.54 


LD 2 


2.24 


3.68 


0.193 


3.19 


1.34 


0.274 


36.94 


1.41 


3.80 


LD 2 


2.24 


3.68 


0.221 


3.15 


1.33 


0.274 


58.50 


2.38 


7.07 


LD 2 


2.24 


3.68 


0.248 


3.12 


1.32 


0.274 


40.19 


1.56 


4.88 


LD 2 


2.24 


3.68 


0.293 


3.06 


1.29 


0.274 


8.69 


0.31 


1.05 


LD 2 


2.14 


3.17 


0.314 


2.17 


0.99 


0.630 


24.32 


0.98 


1.78 


LD 2 


2.14 


3.17 


0.346 


2.14 


0.98 


0.630 


60.74 


2.25 


4.07 


LD 2 


2.14 


3.17 


0.367 


2.11 


0.97 


0.630 


83.29 


3.10 


5.98 


LD 2 


2.14 


3.17 


0.385 


2.09 


0.96 


0.630 


89.87 


3.22 


7.08 


LD 2 


2.14 


3.17 


0.404 


2.07 


0.95 


0.630 


66.90 


2.65 


5.57 


LD 2 


2.14 


3.17 


0.424 


2.05 


0.94 


0.630 


48.89 


1.85 


4.44 


LD 2 


2.14 


3.17 


0.451 


2.02 


0.92 


0.630 


27.91 


1.18 


2.74 


LD 2 


2.14 


3.17 


0.481 


1.98 


0.91 


0.630 


18.00 


0.93 


1.99 


LD 2 


4.74 


4.06 


0.507 


4.48 


1.55 


0.264 


2.76 


0.10 


0.34 


LD 2 


4.74 


4.06 


0.592 


4.40 


1.52 


0.264 


5.54 


0.22 


0.67 


LD 2 


4.74 


4.06 


0.668 


4.34 


1.50 


0.264 


5.31 


0.22 


0.65 


LD 2 


3.94 


3.98 


0.359 


4.18 


1.52 


0.391 


2.47 


0.09 


0.27 


LD 2 


3.94 


3.98 


0.393 


4.15 


1.51 


0.391 


6.20 


0.23 


0.57 


LD 2 


3.94 


3.98 


0.417 


4.12 


1.50 


0.391 


9.38 


0.35 


0.77 


LD 2 


3.94 


3.98 


0.442 


4.10 


1.49 


0.391 


11.94 


0.44 


0.95 


LD 2 


3.94 


3.98 


0.469 


4.08 


1.48 


0.391 


10.97 


0.42 


0.89 


LD 2 


3.94 


3.98 


0.496 


4.05 


1.48 


0.391 


10.04 


0.41 


0.85 


LD 2 


3.94 


3.98 


0.525 


4.03 


1.47 


0.391 


8.22 


0.31 


0.76 



TABLE VII: Extracted cross sections and their uncertainties for deuterium target at eight kinematics settings. The Q 2 , W, 
— t, Pit and e values for each kinematics are shown at the same time. 
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Target 


Q 


W 


-t 




pUM 


e 


da 
dtdr >( -' AJ 


stat. err. 


sys. err. 




GeV 2 


GeV 


GeV^ 


GeV/c 


GeV/c 






lib/GeV s 


fib/GeV' 6 


C 


2.24 


14.01 


0.159 


3.23 


2.51 


0.558 


26.45 


1.02 


2.02 


c 


2.24 


14.01 


0.177 


3.21 


2.49 


0.558 


41.07 


1.61 


2.94 


c 


2.24 


14.01 


0.189 


3.20 


2.48 


0.558 


50.06 


1.76 


3.26 


c 


2.24 


14.01 


0.200 


3.18 


2.47 


0.558 


54.14 


2.09 


3.35 


c 


2.24 


14.01 


0.209 


3.17 


2.46 


0.558 


49.20 


1.88 


3.28 


c 


2.24 


14.01 


0.218 


3.16 


2.45 


0.558 


57.30 


2.12 


3.51 


c 


2.24 


14.01 


0.227 


3.14 


2.44 


0.558 


53.93 


2.01 


3.37 


c 


2.24 


14.01 


0.237 


3.13 


2.43 


0.558 


56.09 


2.08 


3.35 


c 


2.24 


14.01 


0.246 


3.12 


2.42 


0.558 


51.33 


1.92 


3.17 


c 


2.24 


14.01 


0.257 


3.10 


2.41 


0.558 


47.62 


1.77 


2.95 


c 


2.24 


14.01 


0.271 


3.08 


2.39 


0.558 


42.88 


1.48 


2.84 


c 


4.04 


14.59 


0.388 


3.96 


2.91 


0.257 


4.79 


0.32 


0.61 


c 


4.04 


14.59 


0.498 


3.85 


2.83 


0.257 


7.08 


0.36 


0.91 


c 


3.02 


14.22 


0.261 


3.49 


2.65 


0.450 


14.12 


0.55 


1.77 


c 


3.02 


14.22 


0.315 


3.43 


2.61 


0.450 


20.13 


0.66 


2.57 


c 


3.02 


14.22 


0.367 


3.37 


2.56 


0.450 


22.06 


0.69 


2.78 


c 


1.12 


13.67 


0.071 


2.78 


2.24 


0.504 


95.54 


3.80 


9.86 


c 


1.12 


13.67 


0.079 


2.77 


2.23 


0.504 


231.96 


8.73 


25.10 


c 


1.12 


13.67 


0.085 


2.76 


2.22 


0.504 


253.64 


9.85 


29.19 


c 


1.12 


13.67 


0.089 


2.75 


2.21 


0.504 


264.99 


9.82 


31.33 


c 


1.12 


13.67 


0.094 


2.74 


2.20 


0.504 


256.93 


10.00 


30.78 


c 


1.12 


13.67 


0.099 


2.73 


2.19 


0.504 


248.65 


9.55 


30.26 


c 


1.12 


13.67 


0.104 


2.72 


2.19 


0.504 


229.94 


9.41 


26.27 


c 


1.12 


13.67 


0.109 


2.71 


2.18 


0.504 


184.15 


8.36 


20.85 


c 


1.12 


13.67 


0.113 


2.70 


2.17 


0.504 


172.23 


8.24 


16.81 


c 


1.12 


13.67 


0.119 


2.68 


2.16 


0.504 


145.60 


6.24 


14.24 


c 


1.12 


13.67 


0.125 


2.67 


2.15 


0.504 


129.83 


6.11 


10.75 


c 


1.12 


13.67 


0.133 


2.65 


2.13 


0.504 


99.01 


4.60 


7.17 


c 


2.24 


14.01 


0.181 


3.21 


2.49 


0.275 


37.50 


1.60 


3.05 


c 


2.24 


14.01 


0.219 


3.16 


2.45 


0.275 


44.01 


1.85 


3.75 


c 


2.24 


14.01 


0.251 


3.11 


2.41 


0.275 


45.06 


1.84 


3.74 


c 


2.20 


13.17 


0.292 


2.17 


1.77 


0.642 


27.97 


1.35 


3.47 


c 


2.20 


13.17 


0.335 


2.13 


1.73 


0.642 


45.35 


1.84 


5.66 


c 


2.20 


13.17 


0.363 


2.09 


1.71 


0.642 


55.11 


2.75 


6.97 


c 


2.20 


13.17 


0.386 


2.07 


1.69 


0.642 


56.99 


2.25 


7.26 


c 


4.73 


15.03 


0.397 


4.59 


3.27 


0.263 


2.42 


0.16 


0.30 


c 


4.73 


15.03 


0.500 


4.50 


3.21 


0.263 


4.04 


0.19 


0.51 


c 


3.94 


14.77 


0.320 


4.23 


3.08 


0.391 


5.49 


0.22 


0.35 


c 


3.94 


14.77 


0.378 


4.17 


3.04 


0.391 


8.70 


0.34 


0.55 


c 


3.94 


14.77 


0.429 


4.12 


3.00 


0.391 


10.04 


0.34 


0.62 



TABLE VIII: Extracted cross sections and their uncertainties for carbon target at eight kinematics settings. The Q 2 , W, — t, 
and e values for each kinematics are shown at the same time. 
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Target 


Q 


W 


-t 




pUM 

IT 


e 


da 
dtdr >( -' AJ 


stat. err. 


sys. err. 




GeV" 


GeV 


GeV" 


GeV/c 


GeV/c 






lib/GeV s 


fib/GeV' 6 


Cu 


2.24 


62.36 


0.161 


3.23 


3.05 


0.558 


81.49 


3.06 


6.15 


Cu 


2.24 


62.36 


0.183 


3.21 


3.03 


0.558 


112.51 


4.20 


8.81 


Cu 


2.24 


62.36 


0.197 


3.19 


3.01 


0.558 


118.28 


4.48 


9.30 


Cu 


2.24 


62.36 


0.208 


3.17 


2.99 


0.558 


131.91 


5.07 


9.97 


Cu 


2.24 


62.36 


0.219 


3.16 


2.98 


0.558 


131.81 


4.85 


9.93 


Cu 


2.24 


62.36 


0.230 


3.14 


2.96 


0.558 


134.78 


4.92 


9.81 


Cu 


2.24 


62.36 


0.241 


3.12 


2.95 


0.558 


125.55 


4.81 


9.23 


Cu 


2.24 


62.36 


0.253 


3.11 


2.93 


0.558 


122.13 


4.47 


8.79 


Cu 


2.24 


62.36 


0.265 


3.09 


2.91 


0.558 


112.68 


4.25 


8.12 


Cu 


2.24 


62.36 


0.278 


3.07 


2.90 


0.558 


106.16 


4.28 


7.61 


Cu 


4.04 


63.12 


0.364 


3.98 


3.71 


0.257 


10.98 


0.86 


1.51 


Cu 


4.04 


63.12 


0.465 


3.89 


3.62 


0.257 


17.60 


0.81 


2.37 


Cu 


3.02 


62.64 


0.250 


3.51 


3.29 


0.450 


32.07 


1.37 


3.98 


Cu 


3.02 


62.64 


0.305 


3.45 


3.24 


0.450 


46.62 


1.40 


5.01 


Cu 


3.02 


62.64 


0.363 


3.38 


3.17 


0.450 


47.83 


1.61 


4.74 


Cu 


1.12 


61.93 


0.070 


2.79 


2.66 


0.503 


122.69 


4.62 


9.51 


Cu 


1.12 


61.93 


0.083 


2.76 


2.63 


0.503 


618.86 


18.04 


49.79 


Cu 


1.12 


61.93 


0.091 


2.75 


2.62 


0.503 


614.14 


19.83 


48.08 


Cu 


1.12 


61.93 


0.098 


2.73 


2.60 


0.503 


552.29 


18.40 


43.14 


Cu 


1.12 


61.93 


0.105 


2.71 


2.58 


0.503 


542.30 


16.21 


39.69 


Cu 


1.12 


61.93 


0.115 


2.69 


2.56 


0.503 


484.09 


14.21 


33.64 


Cu 


1.12 


61.93 


0.126 


2.67 


2.54 


0.503 


380.04 


12.93 


25.42 


Cu 


2.24 


62.35 


0.186 


3.19 


3.01 


0.282 


109.18 


4.45 


9.03 


Cu 


2.24 


62.35 


0.241 


3.11 


2.94 


0.282 


109.87 


4.09 


8.17 


Cu 


4.73 


63.70 


0.468 


4.53 


4.18 


0.264 


10.76 


0.38 


1.40 


Cu 


4.73 


63.70 


0.647 


4.36 


4.02 


0.264 


10.77 


0.64 


1.42 


Cu 


3.94 


63.34 


0.314 


4.24 


3.93 


0.391 


14.16 


0.62 


1.77 


Cu 


3.94 


63.34 


0.385 


4.17 


3.87 


0.391 


21.61 


0.74 


2.00 


Cu 


3.94 


63.34 


0.450 


4.10 


3.80 


0.391 


23.31 


0.75 


1.94 



TABLE IX: Extracted cross sections and their uncertainties for copper target at seven kinematics settings. The Q 2 , W, — t, Pn 
and e values for each kinematics are shown at the same time. 



Target 


Q 


W 


-t 


Pit 


pCM 


e 


da 


stat. err. 


sys. err. 




GeV 2 


GeV 


GeV" 


GeV/c 


GeV/c 




lib/GeV' 6 


fib/GeV* 


[ibjCeV 6 


Au 


2.16 


186.68 


0.164 


3.23 


3.17 


0.557 


201.07 


7.93 


24.37 


Au 


2.16 


186.68 


0.221 


3.15 


3.09 


0.557 


247.23 


6.02 


29.46 


Au 


4.03 


187.48 


0.438 


3.92 


3.82 


0.256 


33.28 


3.04 


3.15 


Au 


3.02 


186.97 


0.254 


3.50 


3.43 


0.449 


63.50 


3.79 


9.00 


Au 


3.02 


186.97 


0.323 


3.43 


3.36 


0.449 


90.51 


3.56 


12.80 


Au 


1.12 


186.22 


0.070 


2.79 


2.74 


0.504 


149.26 


6.50 


14.36 


Au 


1.12 


186.22 


0.083 


2.76 


2.72 


0.504 


1052.70 


37.33 


104.4 


Au 


1.12 


186.22 


0.092 


2.74 


2.70 


0.504 


1085.12 


37.20 


109.3 


Au 


1.12 


186.22 


0.101 


2.72 


2.68 


0.504 


1033.13 


36.08 


103.5 


Au 


1.12 


186.22 


0.110 


2.70 


2.66 


0.504 


908.06 


34.59 


85.81 


Au 


1.12 


186.22 


0.118 


2.68 


2.64 


0.504 


881.15 


34.71 


75.30 


Au 


4.73 


188.08 


0.545 


4.46 


4.34 


0.263 


20.08 


0.97 


2.79 


Au 


3.95 


187.70 


0.340 


4.21 


4.11 


0.391 


30.82 


1.38 


4.31 


Au 


3.95 


187.70 


0.424 


4.13 


4.03 


0.391 


42.07 


1.77 


5.68 



TABLE X: Extracted cross sections and their uncertainties for gold target at six kinematics settings. The Q 2 , W, —t, P^ and 
e values for each kinematics are shown at the same time. 
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Target 


Q 


W 


-t 




pUM 


e 


da 
dtdr >( -' AJ 


stat. err. 


sys. err. 




GeV 2 


GeV 


GeV^ 


GeV/c 


GeV/c 




ub/GeV 6 


ub/GeV 6 


ub/GeV' 6 


Al 


2.11 


28.20 


0.172 


3.22 


2.84 


0.554 


100.26 


4.31 


11.57 


Al 


2.11 


28.20 


0.228 


3.14 


2.77 


0.554 


82.38 


3.56 


7.81 


Al 


3.99 


28.87 


0.404 


3.94 


3.37 


0.256 


12.90 


3.32 


1.58 


Al 


2.95 


28.44 


0.310 


3.44 


2.99 


0.448 


38.98 


2.09 


5.05 


Al 


1.11 


27.79 


0.077 


2.78 


2.49 


0.502 


401.69 


15.79 


53.74 


Al 


1.11 


27.79 


0.105 


2.72 


2.44 


0.502 


312.93 


11.48 


27.55 


Al 


2.07 


28.21 


0.186 


3.21 


2.82 


0.267 


82.96 


10.66 


10.93 


Al 


2.09 


27.25 


0.286 


2.15 


1.95 


0.648 


81.05 


9.11 


10.05 


Al 


4.57 


29.41 


0.458 


4.54 


3.81 


0.261 


9.44 


1.86 


1.18 


Al 


3.83 


29.09 


0.384 


4.17 


3.54 


0.389 


15.80 


1.36 


2.04 



TABLE XI: Extracted cross sections and their uncertainties for aluminum target at eight kinematics settings. The Q 2 , W, — t, 
and e values for each kinematics are shown at the same time. 

Appendix C 



TVi root 


O 2 


p 


k 




T 


stat. 


ays. 




stat. 


a. 




stat. 






Vjrc V 


VJC V j c 


Vjrc V j C 
























1.1 


2.8 


0.23 


0.50 


0.98 


0.02 


0.03 
















2.2 


3.2 


0.41 


0.56 


1.01 


0.02 


0.03 














±J LJ2 


o.u 


"\ 4 


U.UU 


n 4^ 




n D9 


n 04 














±J LJ2 




4.1 


U. 1 u 


u.oy 




U.Ui 


n n4 














±j lj 2 


4.7 


4.4 


0.79 


0.26 


1.03 


0.03 


0.04 
















2.2 


3.2 


0.42 


0.27 


1.04 


0.03 


0.03 














r 

J-J.LJ2 


a n 




n 71 

U. 1 ± 




i 07 




n D4 














±J LJ2 


2.2 


2.1 


U.UtJ 


u.uo 




n ri9 


u.uo 














c 


1.1 


2.8 


0.23 


0.50 


0.67 


0.01 


0.02 


0.68 


0.01 


0.02 








C 


2.2 


3.2 


0.41 


0.56 


0.65 


0.01 


0.02 


0.65 


0.01 


0.02 








C 


3.0 


3.4 


0.56 


0.45 


0.68 


0.02 


0.03 


0.68 


0.02 


0.03 








C 


3.9 


4.1 


0.70 


0.39 


0.77 


0.02 


0.03 


0.73 


0.02 


0.03 








C 


4.7 


4.4 


0.79 


0.26 


0.70 


0.03 


0.03 


0.69 


0.03 


0.03 








C 


2.2 


3.2 


0.42 


0.27 


0.60 


0.02 


0.02 


0.57 


0.02 


0.02 








C 


4.0 


3.9 


0.71 


0.25 


0.68 


0.03 


0.03 


0.64 


0.03 


0.03 








C 


2.2 


2.1 


0.65 


0.63 


0.66 


0.02 


0.02 


0.66 


0.02 


0.02 








Al 


1.1 


2.8 


0.23 


0.50 


0.49 


0.01 


0.02 


0.50 


0.02 


0.02 


0.73 


0.02 


0.03 


Al 


2.2 


3.2 


0.41 


0.56 


0.52 


0.02 


0.02 


0.52 


0.02 


0.02 


0.80 


0.03 


0.03 


Al 


3.0 


3.4 


0.56 


0.45 


0.57 


0.03 


0.02 


0.57 


0.03 


0.02 


0.83 


0.05 


0.03 


Al 


3.9 


4.1 


0.70 


0.39 


0.59 


0.05 


0.03 


0.56 


0.05 


0.03 


0.77 


0.07 


0.03 


Al 


4.7 


4.4 


0.79 


0.26 


0.71 


0.14 


0.03 


0.69 


0.14 


0.03 


1.01 


0.20 


0.04 


Al 


2.2 


3.2 


0.42 


0.27 


0.46 


0.06 


0.02 


0.45 


0.06 


0.02 


0.78 


0.10 


0.03 


Al 


4.0 


3.9 


0.71 


0.25 


0.70 


0.16 


0.03 


0.66 


0.15 


0.02 


1.03 


0.24 


0.03 


Al 


2.2 


2.1 


0.65 


0.63 


0.59 


0.07 


0.02 


0.60 


0.07 


0.02 


0.90 


0.10 


0.03 


Cu 


1.1 


2.8 


0.23 


0.50 


0.45 


0.01 


0.01 


0.46 


0.01 


0.02 


0.67 


0.01 


0.02 


Cu 


2.2 


3.2 


0.41 


0.56 


0.45 


0.01 


0.02 


0.45 


0.01 


0.02 


0.70 


0.01 


0.03 


Cu 


3.0 


3.4 


0.56 


0.45 


0.43 


0.01 


0.02 


0.44 


0.01 


0.02 


0.64 


0.02 


0.03 


Cu 


3.9 


4.1 


0.70 


0.39 


0.52 


0.01 


0.02 


0.49 


0.01 


0.02 


0.67 


0.02 


0.03 


Cu 


4.7 


4.4 


0.79 


0.26 


0.53 


0.02 


0.02 


0.51 


0.02 


0.02 


0.75 


0.04 


0.03 


Cu 


2.2 


3.2 


0.42 


0.27 


0.43 


0.02 


0.02 


0.42 


0.01 


0.02 


0.73 


0.03 


0.03 


Cu 


4.0 


3.9 


0.71 


0.25 


0.51 


0.02 


0.02 


0.48 


0.02 


0.02 


0.75 


0.04 


0.03 


Au 


1.1 


2.8 


0.23 


0.50 


0.28 


0.01 


0.01 


0.28 


0.01 


0.01 


0.41 


0.01 


0.02 


Au 


2.2 


3.2 


0.41 


0.56 


0.29 


0.01 


0.01 


0.28 


0.01 


0.01 


0.44 


0.01 


0.02 


Au 


3.0 


3.4 


0.56 


0.45 


0.29 


0.01 


0.01 


0.30 


0.01 


0.01 


0.43 


0.02 


0.02 


Au 


3.9 


4.1 


0.70 


0.39 


0.34 


0.01 


0.02 


0.32 


0.01 


0.01 


0.44 


0.02 


0.02 


Au 


4.7 


4.4 


0.79 


0.26 


0.33 


0.02 


0.02 


0.32 


0.02 


0.01 


0.46 


0.03 


0.02 


Au 


4.0 


3.9 


0.71 


0.25 


0.31 


0.03 


0.02 


0.29 


0.03 


0.01 


0.46 


0.05 


0.02 



TABLE XII: Extracted nuclear transparency and their uncertainties. The Q 2 , P^, k n and e values for each kinematics are 
shown at the same time. Here the T, Ta,i and Ta,i2 are the nuclear transparency formed with hydrogen, deuterium and carbon 
targets, respectively. The Q 2 dependent model uncertainty is 7.6%, 5.7%, 3.5%, 3.8%, and 3.8% for Q 2 = 1.1, 2.1, 3.0, 3.9, 4.7 
(GeV/c) 2 . respectively. 



